TECHNICAL  REPORT 


Defense  Threat  Reduction  Agency 
8725  John  J.  Kingman  Road,  MS 
6201  Fort  Belvoir,  VA  22060-6201 


DTRA-TR-16-77 

Computational  Characterization  of 
Impact  Induced  Multi-Scale 
Dissipation  in  Reactive  Solid 
Composites 


Distribution  Statement  A.  Approved  for  public  release;  distribution  is 
unlimited. 


July  2016 


HDTRA1-10-1-0018 
Keith  A.  Gonthier 


Prepared  by: 

Louisiana  State  University 
Baton  Rouge,  LA  70803 


DESTRUCTION  NOTICE: 


Destroy  this  report  when  it  is  no  longer  needed. 

Do  not  return  to  sender. 

PLEASE  NOTIFY  THE  DEFENSE  THREAT  REDUCTION 
AGENCY,  ATTN:  DTRIAC/  J9STT,  8725  JOHN  J.  KINGMAN  ROAD, 
MS-6201 ,  FT  BELVOIR,  VA  22060-6201 ,  IF  YOUR  ADDRESS 
IS  INCORRECT,  IF  YOU  WISH  IT  DELETED  FROM  THE 
DISTRIBUTION  LIST,  OR  IF  THE  ADDRESSEE  IS  NO 
LONGER  EMPLOYED  BY  YOUR  ORGANIZATION. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data 
sources,  gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other 
aspect  of  this  collection  of  information,  including  suggestions  for  reducing  the  burden,  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information 
Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other 
provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently  valid  OMB  control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1.  REPORT  DATE  (DD-MM-YYYY) 

00-07-2016 


2.  REPORT  TYPE 

Final  Report 


3.  DATES  COVERED  (From  -  To) 

March  2010-August  2015 


4.  TITLE  AND  SUBTITLE 

Computational  Characterization  of  Impact  Induced  Multi-Scale  Dissipation  in 
Reactive  Solid  Composites 


5a.  CONTRACT  NUMBER 

HDTRA1-10-1-0018 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


6.  AUTHOR(S) 

Gonthier,  Keith,  A 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Louisiana  State  University 
Baton  Rouge,  Louisiana,  70803 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 
Defense  Threat  Reduction  Agency 
8725  John  J.  Kingman  Road 
Fort  Belvoir,  VA  22060-6201 


10.  SPONSOR/MONITOR'S  ACRONYM(S) 

DTRA 


11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

DTRA-TR- 16-77 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Distribution  Statement  A.  Approved  for  public  release;  distribution  is  unlimited. 


13.  SUPPLEMENTARY  NOTES 


14.  ABSTRACT 

Mesoscale  and  macroscale  modeling  and  simulation  (M&S)  were  used  to  characterize  shock  induced  dissipation  in  low 
density  granular  explosives  (HMX,  and  HMX-AI  composites).  Emphasis  was  placed  on  characterizing  how  material 
microstructure  (particle  size  and  packing  density)  and  composition  (metal  mass  fraction)  affect  shock  heating  of  the 
explosive  component  which  establishes  shock  sensitivity  and  survivability.  Predicted  variations  in  hot-spot  formation  rate 
with  shock  strength  were  used  to  develop  a  microstructure-dependent  ignition  and  burn  model  for  low  density  HMX.  The 
model  was  used  to  computationally  examine  initiation  mechanisms  resulting  from  low  pressure  shock  input. 


15.  SUBJECT  TERMS 

Low  Density  HMX,  Porous  HMX,  Meso-Scale  Modeling  and  Simulation,  Hot-Spot  Formation,  Shock  Initiation,  Deflagration 
to  Detonation  Transition,  Ignition  and  Burn  Modeling 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

18.  NUMBER 

19a.  NAME  OF  RESPONSIBLE  PERSON 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

ABSTRACT 

OF 

PAGES 

Allen  Dalton 

U 

U 

U 

SAR 

50 

19b.  TELEPHONE  NUMBER  (Include  area  code) 

703-767-4663 

Standard  Form  298  (Rev.  8/98) 
Prescribed  by  ANSI  Sfd.  Z39.18 


UNIT  CONVERSION  TABLE 


U.S.  customary  units  to  and  from  international  units  of  measurement* 


U.S.  Customary  Units 

iviuiupiy  Dy 

International  Units 

'  Divide  by  ' 

Length/  Area/V  olume 

inch  (in) 

2.54 

x  10-2 

meter  (m) 

foot  (ft) 

3.048 

x  10_1 

meter  (m) 

yard  (yd) 

9.144 

x  10_1 

meter  (m) 

mile  (mi,  international) 

1.609  344 

x  103 

meter  (m) 

mile  (nmi,  nautical,  U.S.) 

1.852 

x  103 

meter  (m) 

barn  (b) 

1 

x  10-28 

square  meter  (m2) 

gallon  (gal,  U.S.  liquid) 

3.785  412 

x  10'3 

cubic  meter  (m3) 

cubic  foot  (ft3) 

2.831  685 

x  10-2 

cubic  meter  (m3) 

Mass/Density 
pound  (lb) 

4.535  924 

x  10_1 

kilogram  (kg) 

unified  atomic  mass  unit  (amu) 

1.660  539 

x  10-27 

kilogram  (kg) 

pound-mass  per  cubic  foot  (lb  ft-3) 

1.601  846 

x  101 

kilogram  per  cubic  meter  (kg  m"3) 

pound-force  (lbf  avoirdupois) 

Energy/W  ork/Power 

4.448  222 

x  10-19 

newton  (N) 

electron  volt  (eV) 

1.602  177 

joule  (J) 

erg 

1 

x  10-7 

joule  (J) 

kiloton  (kt)  (TNT  equivalent) 

4.184 

x  1012 

joule  (J) 

British  thermal  unit  (Btu) 

(thermochemical) 

1.054  350 

x  103 

joule  (J) 

foot-pound-force  (ft  lbf) 

1.355  818 

joule  (J) 

calorie  (cal)  (thermochemical) 

Pressure 

4.184 

joule  (J) 

atmosphere  (atm) 

1.013  250 

x  105 

pascal  (Pa) 

pound  force  per  square  inch  (psi) 

Temperature 

6.984  757 

x  103 

pascal  (Pa) 

degree  Fahrenheit  (°F) 

[T(°F)-32]/1.8 

degree  Celsius  (°C) 

degree  Fahrenheit  (°F) 

[T(°F) +  459.67]/!. 8 

kelvin  (K) 

Radiation 

curie  (Ci)  [activity  of  radionuclides] 

3.7 

x  1010 

per  second  (s_1)  [becquerel  (Bq)] 

roentgen  (R)  [air  exposure] 

2.579  760 

x  10~4 

coulomb  per  kilogram  (C  kg~') 

rad  [absorbed  dose] 

1 

x  10  2 

joule  per  kilogram  (J  kg-1)  [gray  (Gy)] 

rem  [equivalent  and  effective  dose] 

1 

x  10  2 

joule  per  kilogram  (J  kg-1)  [sievert  (Sv)] 

Specific  details  regarding  the  implementation  of  SI  units  may  be  viewed  at  http  ://ww w.bipm.org/en/ si/. 

^Multiply  the  U.S.  customary  unit  by  the  factor  to  get  the  international  unit.  Divide  the  international  unit  by  the  factor  to  get  the 
U.S.  customary  unit. 


2015-11-16 


FINAL  REPORT 


Defense  Threat  Reduction  Agency 


Computational  Characterization  of  Impact  Induced 
Multi-Scale  Dissipation  in  Reactive  Solid  Composites 


Principal  Investigator:  Keith  A.  Gonthier 
Mechanical  and  Industrial  Engineering  Department 
Louisiana  State  University 
Baton  Rouge,  Louisiana  70803 
Email:  gonthier@me.lsu.edu 


Award  Number:  HDTRAl-10-1-0018 
Performance  Period:  March  2010-August  2015 


February  24,  2016 


Contents 


1  Abstract  4 

2  Shock  Dissipation  in  Low  Density  Granular  Explosives  5 

2.1  Significance .  5 

2.2  Goals  and  Methods  .  5 

2.3  Work  Accomplished  .  7 

3  Analysis  of  Hot-Spots  in  Low  Density  Granular  Explosives  14 

3.1  Significance .  14 

3.2  Goals  and  Methods .  15 

3.3  Work  Accomplished  .  15 

4  Ignition  and  Burn  Modeling  for  Low  Density  HMX  28 

4.1  Significance .  28 

4.2  Goals  and  Methods .  28 

4.3  Work  Accomplished  .  28 

5  Personnel  Supported  39 

6  Publications  and  Conference  Presentations  40 

7  Workshop  and  Meeting  Presentations  42 

8  Manuscripts  in  Review/Preparation  42 

9  Abstracts  Accepted  for  Presentation  42 


1 


List  of  Figures 


1  Initial  microstructure  for  an  ensemble  of  explosive  and  metal  particles  having  Xm  = 

0.2  and  (j)so  =  0.835.  Explosive  particles  are  colored  red  and  metal  particles  are 
colored  blue .  7 

2  Predicted  variation  in  Hugoniot  curves  with  \m  {cj)s o  =  0.835):  (a)  Up-D  plane;  (b) 

Up-P  plane;  and  (c)  Up-(j)  plane .  8 

3  Predictions  for  effective  component-specific  spatial  profiles  within  a  (a)  strength  dom¬ 
inated  shock  {Up  =  100  m/s)  and  (b)  pressure  dominated  shock  ( Up  =  500  m/s).  The 
metal  mass  fraction  is  Xm  =  0.20.  The  effective  phase-specific  pressure  and  Von- 
Mises  stress  is  shown  in  each  case.  Standard  deviations  about  the  effective  profiles 

are  shown  as  thin  dashed  curves.  Black  curves  correspond  to  neat  HMX .  10 

4  Predictions  for  hot-spot  mass  fraction  curves  within  a  (a)  strength  dominated  shock 
{Up  =  100  m/s)  and  (b)  pressure  dominated  shock  {Up  =  500  m/s).  The  metal  mass 
fraction  is  Xm  =  0.20.  The  black  curves  are  the  effective  dissipative  heating  rates  of 

the  (a)  explosive  and  (b)  metal  components .  12 

5  Illustration  of  the  temperature  thresholding  technique  used  to  identify  hot-spots. 

Here,  Tth  =  500  K  denotes  the  temperature  threshold .  16 

6  Representative  initial  particle  ensembles  used  with  mesoscale  M&S  to  assess  the 

effect  of  packing  density  on  hot-spot  formation  in  low  density  HMX .  17 

7  Predicted  variation  in  (a)  hot-spot  number  density,  (b)  hot-spot  volume  fraction,  and 

(c)  hot-spot  specific  surface  area  for  each  ensemble  with  piston  speed.  The  dashed 
lines  represent  best-fit  exponential  functions .  18 

8  An  illustration  of  hot-spot  axial  spacing  used  to  compute  the  formation  time  between 

successive  hot-spots  based  on  their  centroid  locations .  19 

9  Predicted  variation  in  hot-spot  formation  time  distribution  with  Up  for  steady  waves 

in  porous  HMX:  (a)  </>0  =  0.68;  (b)  (j)0  =  0.84 .  20 

10  Predicted  variation  in  the  effective  hot-spot  formation  rate  with  (a)  piston  speed  and 

(b)  shock  pressure  for  low  density  HMX .  21 

11  (a)  Predicted  variation  in  the  global  ignition  time  with  wave  strength  for  porous 

HMX  having  different  values  of  <f>0.  (b)  Comparison  of  the  predicted  ignition  curve 
with  initiation  time  data  for  0O  =  0.68 .  23 

12  Initial  particle  ensembles  with  varying  packing  densities  (ps  and  composition  Xm . 

Here,  red  particles  represent  Al,  and  yellow  particles  represent  HMX .  24 

13  Predicted  variations  in  effective  shock  end  states  as  function  of  4>s  and  Am:  (a )P-Up\ 

(b)  D-Up,  (c)  Wpe-Up-  and  (d)  Wpe-P .  25 

14  Predicted  variations  in  the  effective  plastic  work  of  the  explosive  component  as  a 

function  of  cj>s,  and  Xm .  26 

15  Meso-structures  with  varying  particle  sizes,  particle  distributions  and  composition. 

Here,  red  particles  represent  metal,  and  yellow  particles  represent  explosive .  26 

16  Predicted  variations  in  explosive  dissipative  work  with  particle  size  ratio  and  shock 

pressure  for  a  granular  HMX-A1  composite  having  cj)s  =  0.84:  (a)Am  =  0.20  and 
(b)Am  =  0.56 .  27 

17  Predicted  variations  in  hot-spot  volumetric  properties  with  particle  size  ratio  and 
shock  pressure  for  a  granular  HMX-A1  composite  having  (j)s  =  0.84  and  Xm  =  0.20: 

(a)  number  density  and  (b)  volume  fraction .  27 


2 


18  A  comparison  between  mesoscale  and  continuum  theory  predictions  for  shock  dis¬ 
sipation  rate  and  dissipative  work  profiles  corresponding  to  Up  =  500  m/s.  Shock 
profiles  are  shown  for  different  values  of  initial  solid  volume  fraction  </>0.  The  shock 

is  propagating  from  left  to  right  in  the  figure .  31 

19  A  comparison  between  mesoscale  and  continuum  theory  predictions  for  (a)  the  vari¬ 

ation  in  effective  dissipative  work  with  shock  pressure  and  (b)  the  variation  in  com¬ 
paction  shock  rise  time  with  supporting  piston  speed .  31 

20  A  comparison  of  mesoscale  M&S  and  continuum  compaction  theory  predictions  for 
dissipative  work  profiles  induced  by  successive  planar  shock  loading  of  low  density 
HMX  having  (j>o  =  0.68.  The  successive  supporting  piston  speeds  are  JJp\  =  300 
m/s  and  Up2  =  500  m/s.  (a)  Mesoscale  and  (b)  continuum  theory  spatial  fields  at 
various  times.  Fields  are  plotted  in  a  right  propagating  piston- attached  reference 
frame;  the  piston  surface  is  located  at  £  =  0  mm.  (c)  Mesoscale  and  (d)  continuum 


theory  predictions  plotted  in  pressure-dissipated  work  phase  space .  32 

21  Plots  showing  the  (a)  ignition  and  initiation  times,  and  (b)  initiation  distances  for 

</> o  =  0.68  and  <f>0  =  0.83 .  36 

22  Predicted  (a)  bulk/mixture  pressure  contours  showing  compaction  shock  and  igni¬ 
tion  front  trajectories  in  the  characteristic  plane;  (b)  solid/explosive  volume  fraction 
contours  in  the  characteristic  plane;  (c)  variation  in  bulk  pressure  with  position  at 
various  times;  and  (d)  variation  in  dissipative  compaction  work  with  position  at 
various  times.  The  material  is  low  density  HMX  (</> o  =  0.68)  and  the  input  shock 


corresponds  to  a  supporting  piston  speed  of  Up  =  200  m/s.  Position  is  shown  in  a 
piston-attached  reference  frame .  38 


3 


Computational  Characterization  of  Impact  Induced 
Multi-Scale  Dissipation  in  Reactive  Solid  Composites 


Keith  A.  Gonthier 
Associate  Professor 

Mechanical  and  Industrial  Engineering  Department 
Louisiana  State  University 
Baton  Rouge,  Louisiana  70803 
Phone:  (225)  578-5915 
FAX:  (225)  578-5924 
Email:  gonthier@me.lsu.edu 


1  Abstract 

Reactive  solid  composites  often  consist  of  mixtures  of  high-explosive  and  metal  particles  (size  ~ 
0.001-200  /-ini),  and  polymeric  binder.  It  remains  fundamentally  unclear,  even  for  idealized  systems, 
how  microstructure  (particle  size,  shape,  and  packing),  component  thermomechanical  properties, 
and  metal  and  binder  mass  fractions,  affect  impact  induced  heating  of  the  explosive  component 
which  establishes  their  impact  sensitivity  and  survivability. 

This  modeling  and  computational  study  had  three  key  objectives:  1)  To  examine  mesoscale  dis¬ 
sipation  by  compaction  shocks  in  low  density  granular /particulate  explosives  (HMX,  and  HMX-A1 
composites).  For  this  objective,  mesoscale  modeling  and  simulation  was  performed  using  a  tech¬ 
nique  that  is  well-suited  for  describing  both  volumetric  and  surface  dissipation  in  granular  systems 
that  can  induce  thermally  activated  phenomena.  The  technique  incorporates  a  thermoelastic- 
viscoplastic  and  stick-slip  friction  theory  for  each  component  to  describe  nonlinear  deformation 
and  motion,  inter-particle  friction,  and  plastic  work.  The  relative  importance  of  volumetric  and 
surface  dissipation  within  shock  profiles  was  characterized,  and  its  dependence  on  shock  strength 
and  metal  mass  fraction  was  examined.  2)  To  characterize  how  the  microstructure  and  composition 
of  low  density  granular  explosives  (HMX,  and  HMX-A1  composites)  affect  shock  induced  formation 
of  hot-spots.  For  this  objective,  inert  temperature  field  predictions  at  the  particle  scale  were  com¬ 
bined  with  a  thresholding  strategy  to  identify  hot-spot  fields  and  to  characterize  their  distributions 
in  intensity,  geometry,  and  spatial  proximity  in  the  deformed  material  configuration  behind  shocks. 
Such  distributions  are  significant  because  they  establish  local  ignition  and  control  the  rate  of  flame 
spread  within  the  material.  3)  To  develop  a  thermodynamically  compatible  macroscale  ignition  and 
burn  model  that  explicitly  incorporates  computationally  derived  relations  between  microstructure, 
shock  strength,  and  hot-spots.  The  model  was  used  to  computationally  examine  how  shock  induced 
transition  to  detonation  in  low  density  HMX  is  affected  by  input  shock  strength  and  initial  packing 
density. 
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This  report  briefly  summarizes  the  significance,  key  goals  and  methods,  and  work  accomplished 
for  each  objective  of  this  study.  Only  brief  discussions  of  key  results  are  given;  detailed  expositions 
can  be  found  in  the  cited  articles. 

2  Shock  Dissipation  in  Low  Density  Granular  Explosives 

2.1  Significance 

Mesoscale  Modeling  and  Simulation  (M&S)  is  a  potentially  useful  tool  to  characterize  shock  induced 
dissipation  within  granular  explosives  that  drives  thermally  activated  phenomena  responsible  for 
ignition  and  flame  spread.  The  utility  of  this  approach  is  in  its  ability  to  resolve  thermomechanical 
fields  at  the  particle  scale  that  are  difficult  to  experimentally  interrogate.  Therefore,  mesoscale 
M&S  can  principally  provide  information  about  the  influence  of  microstructure  on  the  material’s 
shock  behavior.  This  information  can  be  potentially  used  to  engineer  microstructures  for  desired 
performance,  and  may  also  be  used  to  guide  development  of  microstructure-dependent  continuum 
theories  for  shock  induced  ignition  and  burn  of  granular  explosives. 

2.2  Goals  and  Methods 

The  goal  of  this  effort  was  to  computationally  examine  dissipation  within  and  behind  uniaxial  com¬ 
paction  shocks  in  low  density  granular  explosives  (HMX,  and  HMX-Al  composites).  Mesoscale  M&S 
was  performed  to  identify  how  A1  mass  fraction  affects  the  spatial  structure  of  piston  supported 
compaction  shocks.  An  identical  particle  size  distribution  and  packing  arrangement  was  used  for 
all  simulations  to  isolate  the  effects  of  material  composition  (i.e. ,  A1  mass  fraction)  from  those  due 
to  variations  in  microstructure  and  packing  arrangement.  This  study  represents  a  systematic  ex¬ 
tension  of  the  work  by  Panchadhara  and  Gonthier  [11]  which  strictly  focused  on  compaction  shocks 
in  low  density  HMX  to  establish  a  benchmark  for  characterizing  the  effects  of  metal  additives  on 
shock  behavior.  Emphasis  was  placed  in  this  effort  on  examining  how  the  inclusion  of  A1  affects 
explosive  component  heating  and  hot-spots  relative  to  that  of  neat  HMX.  Here,  neat  material  refers 
to  one  that  consists  of  only  a  single  component. 

The  physical  problem  was  posed  as  a  multi-particle  contact  problem  for  the  evolution  of  ther- 
momechanical  fields  within  particles.  The  coupled  initial  boundary  value  problem  (IBVP)  for  the 
displacement  field  u  and  temperature  field  T  within  explosive  and  metal  particles  was  locally  de¬ 
scribed  by  evolution  equations  for  linear  momentum  and  thermal  energy.  All  particles  were  initially 
stationary  and  unstressed.  A  hyperelastic,  multiplicative,  finite  strain  constitutive  theory  was  used 
to  model  the  stress-strain  behavior  of  explosive  and  metal  particles  resulting  from  impact.  A 
Perzyna  type  over-stress  model,  coupled  with  an  associative  flow  rule  and  a  Von-Mises  type  yield 
criterion,  was  used  to  describe  the  evolution  of  inelastic  strain.  The  hyperelastic  material  behavior 
was  formulated  in  terms  of  the  Kirchhoff  stress  tensor  r  =  Ya=i  (Ta)  n*  ®  n?:>  where  n*  and  Ta  are 
its  orthonormal  eigenvectors  and  corresponding  eigenvalues  which  were  computed  from  a  stored 
strain  energy  function  'I'  by  Ta  =  <9T/<91n(Ae,i).  Here,  Ae,i  are  the  eigenvalues  of  the  elastic  com¬ 
ponent  of  the  deformation  gradient  F  =  1  +  Vou,  where  Vo  =  <9(-)/3X  is  the  Lagrangian  gradient 
operator.  The  Kirchhoff  stress  tensor  is  related  to  the  Cauchy  stress  tensor  by  r  =  Jer,  which  can 
be  decomposed  into  volumetric  and  deviatoric,  components:  Ta  =  J  p  +  fa,  where  the  pressure  is 
given  by  p  =  d^/dln (J),  the  stress  deviator  is  given  by  fa  =  3T/31n(A e,i)i  and  A e,i  = 

The  Jacobian  J  =  det(F)  characterizes  volumetric  deformation. 
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The  contact  traction  vector  between  particles  was  prescribed  by  tc  =  tjyn  +  tt,  where  fjy  is  the 
local  compressive  normal  traction  on  contact  boundaries  defined  by  =  —  n  cr  n  and  n  is  the  unit 
outward  normal  vector.  The  frictional  traction  vector  tj  was  estimated  using  Amontons-Coulomb 
law  of  dry  friction  for  which  particle  stick  or  slip  was  possible  depending  on  the  system  dynamics. 
A  contact  was  assumed  to  stick  if  fa  =  |t*|  —  ptjy  5s  0,  where  p,  is  the  limiting  friction  coefficient; 
if  the  contact  was  sliding,  then  the  tangential  traction  was  given  by  |tt|  =  The  local  heat 

flux  due  to  frictional  heating  was  given  by  qj  =  tc  ■  vr,  where  vr  is  the  relative  tangential  velocity. 

The  total  heat  flux  was  distributed  between  the  contacting  bodies  according  to  qj  1  =  uqj  and 
qj  =  (1  —  cu)qf,  where  the  partition  function  uj  was  given  by 


The  notation  superscript  (1)  and  (2)  represent  properties  associated  with  the  two  particles  in 
contact  which  consisted  of  either  two  explosive  particles,  a  metal  particle  and  an  explosive  particle, 
or  two  metal  particles.  Thermal  contact  resistance  between  particles  was  ignored. 

An  additive  form  of  the  stored  energy  function  was  used  to  describe  the  behavior  of  HMX  and 
Al: 

*  =  [ln( Ae^ )] 2  +  [ln( J)]2  -  ^(T  -  T0)ln(J)  +  p0cv  [(T  -T0)-T  In  .  (2) 

Here,  G  is  the  shear  modulus,  k  is  the  bulk  modulus,  a  is  the  linear  coefficient  of  thermal  expansion, 
and  To  is  the  ambient  temperature  of  the  initially  stress  free  ensemble.  The  effective  plastic  strain 
rate  was  estimated  using  a  Perzyna  type  overstress  model  ep  =<  d/ 7  >,  where  the  overstress 
d  =  Te  —  To,  7  is  the  plastic  viscosity  parameter,  and  <  •  >  are  the  Macaulay  brackets.  In  this 
expression  re  =  \J\t  :  r  is  the  effective  (or  Von  Mises)  stress  and  to  is  the  initial  yield  stress. 
Explicit  relations  for  pressure  and  deviatoric  stress  are  given  by 

P  =  j  ^fcln(J)  -  3ka(T  -  T0)— — ^  ,  fa  =  2Gln(Aeji).  (3) 

Key  differences  exist  in  several  material  properties  for  HMX  and  Al  that  influenced  the  defor¬ 
mation  response  of  the  mixture.  First,  the  characteristic  acoustic  impedance,  given  by  Z  =  poCi, 
where  po  and  Ci  are  the  ambient  density  and  longitudinal  sound  speed,  is  appreciably  higher  for  the 
metal  than  the  explosive  (. Zmetai/Zexpi  ps  2.3).  Because  the  acoustic  pressure  goes  as  P  ~  Zu.  where 
u  is  the  longitudinal  particle  speed,  the  inclusion  of  metal  increased  the  average  (or  effective)  wave 
pressure  over  that  of  neat  granular  HMX  for  a  given  microstructure  and  impact  speed.  An  increase 
in  effective  wave  pressure  generally  facilitated  the  elimination  of  porosity  and  enhanced  hot-spot 
formation  within  the  explosive  component,  though  the  degree  of  enhancement  was  mitigated  by 
the  extent  of  inelastic  metal  deformation  because  of  its  low  yield  stress  (To)me4a; / TQjexp;  ps  0.67). 
Second,  the  thermal  conductivity  of  the  metal  is  significantly  higher  than  that  of  the  explosive 
(kr, metal / kr,expi  ~  474)  which  affected  the  partitioning  of  frictionally  dissipated  energy  at  metal- 
explosive  contact  surfaces.  Based  on  ambient  material  densities,  Eq.  (1)  gives  uj  «  0.05;  thus, 
approximately  95%  of  frictionally  dissipated  energy  at  contacts  was  absorbed  by  the  metal  which, 
in  part,  suppressed  explosive  hot-spot  formation.  The  extent  of  suppression  was  dependent  on  both 
metal  mass  fraction  and  microstructure. 
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Figure  1:  Initial  microstructure  for  an  ensemble  of  explosive  and  metal  particles  having  Am  =  0.2 
and  (j)so  =  0.835.  Explosive  particles  are  colored  red  and  metal  particles  are  colored  blue. 


The  computational  technique  used  the  finite-element  method  (FEM),  coupled  with  a  radial 
return  stress  update  algorithm,  to  numerically  integrate  the  time-dependent,  2-D  conservation 
principles  and  viscoplastic  flow  rule  governing  deformation  of  individual  particles,  and  used  the 
discrete-element  method  (DEM)  to  account  for  interactions  between  particles.  The  DEM  used 
a  distributed,  conservative  potential  based  penalty  method  whereby  the  normal  contact  traction 
between  particles  was  estimated  by  penalizing  their  penetration,  and  frictional  tractions  were  esti¬ 
mated  by  a  penalty  regularized  Amontons  Coulomb  law.  Particles  were  discretized  using  constant 
strain,  triangular  (CST)  finite  elements.  A  temporally  second-order  accurate,  explicit  numerical 
technique  was  used  to  integrate  the  finite  element  equations  for  nodal  displacements  and  temper¬ 
atures.  The  constitutive  theory  was  separately  validated  using  Hugoniot  data  for  granular  HMX 
and  porous  Al. 

2.3  Work  Accomplished 

Shock  structures  were  computationally  examined  based  on  effective  shock  pressure  relative  to  the 
material  crush-up  strength.  Work  accomplished  as  part  of  this  effort  along  with  key  observations 
are  briefly  summarized  below.  Publications  and  presentations  resulting  from  this  work  are  indicated 
by  the  cited  references. 

1.  Analysis  of  computed  Hugoniots.  [Chakravarthy,  et.  ah,  2010;  Chakravarthy  and  Gonthier, 
2010;  Gonthier,  2010;  Chakravarthy  and  Gonthier,  2011;  Chakravarthy,  et  al.,  2013  A; 
Chakravarthy,  2014],  The  deformation  response  of  material  having  an  initial  solid  volume 
fraction  of  =  0.835  (porosity  1  —  (f>®  =  0.165)  was  characterized  for  different  shock  pres¬ 
sures  and  Al  mass  fractions.  The  microstructure  consisted  of  initially  circular,  randomly 
packed  HMX  and  Al  particles  having  three  discrete  particle  sizes  of  d  =  (40,  50,  60)  /irn,  with 
a  mean  of  d  =  50  gm,  as  illustrated  in  Fig.  1  for  a  particle  ensemble  having  an  Al  mass 
fraction  of  Am  =  0.2. 

Figure  2  summarizes  predicted  Hugoniot  curves  for  Xm  =  0,  0.20,  0.56,  and  1  in  the  Up-D  and 
Up-P  planes,  respectively,  where  P  =  ( Pe(f>e  +  Pm4>m)/4’s  is  the  effective  equilibrium  pressure 
behind  the  shock,  (f>s  =  cf>e  +  </>m  is  the  total  solid  volume  fraction  behind  the  shock,  and  D 
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Figure  2:  Predicted  variation  in  Hugoniot  curves  with  Am  (< f>so  =  0.835):  (a)  Up-D  plane;  (b)  Up-P 
plane;  and  (c)  Up-(f>  plane. 


is  the  corresponding  shock  speed.  Here,  piston  speeds  within  the  range  50  <  Up  <  500  m/s 
were  considered,  where  speeds  of  approximately  Up  >  250  m/s  resulted  in  nearly  complete 
material  crush-up  and  consolidation  (cj)s  ~  1)  for  all  values  of  Xm.  This  critical  value  of 
piston  speed  separates  two  well-defined  wave  categories:  a  low  pressure  category  referred  to 
as  strength  dominated  because  material  strength  is  important  in  preventing  complete  pore 
collapse;  and  a  high  pressure  category  referred  to  as  pressure  dominated  because  the  resulting 
pressure  is  sufficient  to  completely  eliminate  porosity.  Material  strength  plays  a  less  significant 
role  than  pressure  for  pressure  dominated  shocks. 

A  few  observations  are  noteworthy.  First,  the  effective  shock  pressure  monotonically  increases 
with  Xm  for  fixed  Up  due  to  the  high  metal  acoustic  impedance  which  causes  reflection  of  acous¬ 
tic  energy  at  explosive-metal  interfaces  thereby  enhancing  the  deformation  and  pressure  of 
explosive  particles.  This  observation  agrees  with  data  for  both  neat  granular  HMX  ( \m  =  0) 
and  neat  porous  aluminum  ( Xm  =  1)  [9].  Moreover,  the  inclusion  of  metal  more  significantly 
affects  shock  pressure  for  pressure  dominated  waves  due  to  enhanced  compression  of  consol¬ 
idated  material.  Second,  variations  in  D  with  Xm  qualitatively  differ  between  the  two  wave 
categories.  For  strength  dominated  shocks,  D  slightly  decreases  with  increasing  Xm  whereas 
the  opposite  trend  is  predicted  for  pressure  dominated  shocks.  Though  speculative,  this  pre¬ 
diction  is  plausible  because  the  lower  metal  yield  strength  can  enhance  dissipation  by  plastic 
work  during  pore  collapse,  even  for  neat  granular  metal  (Am  =  1),  resulting  in  slightly  lower 
wave  speeds  than  those  associated  with  neat  granular  HMX.  However,  dissipation  by  plastic 
flow  becomes  more  geometrically  constrained  within  shock  profiles  for  pressure  dominated 
shocks  due  to  the  elimination  of  porosity  resulting  in  higher  wave  speeds.  The  weak  variation 
in  D  with  Up  for  strength  dominated  shocks  is  largely  due  to  inter-particle  friction  during 
pore  collapse  which  enhances  material  rigidity.  Simulations  performed  for  frictionless  contact 
indicate  a  larger  decrease  in  D  as  Up  — >  0  for  strength  dominated  waves  (as  illustrated  Item 
5  of  this  section).  It  is  possible  that  particle  fracture  and  rearrangement,  and  3D  effects, 
may  cause  a  larger  variation  in  D  with  Up  for  strength  dominated  shocks  by  enhancing  the 
mobility  of  fractured  particles  during  crush- up. 

Figure  2(c)  gives  the  predicted  variation  in  cj)s  and  <j)e  behind  the  shock  with  Up.  These 
predictions  illustrate  that  pressure  dominated  shocks  result  in  nearly  complete  material  con¬ 
solidation  for  all  values  of  Xm.  Moreover,  mild  variations  in  (j>e  with  Up  are  predicted  to 
occur  following  consolidation  (with  </>m  =  1  —  </>e) ,  particularly  for  high  pressure  shocks,  due 
to  differences  in  component-specific  compressibility.  Because  of  the  high  metal  bulk  modulus 
(Kmetai/ Kexpi  «  5.8),  the  explosive  is  more  easily  compressed  under  high  pressure  confine¬ 
ment  resulting  in  a  slight  reduction  in  its  volume  fraction.  This  effect  is  more  pronounced  for 
larger  values  of  Xm  as  shock  pressure  increases. 

2.  The  spatial  structure  of  compaction  shocks  was  computationally  examined.  [Chakravarthy, 
et.  al.,  2010;  Chakravarthy  and  Gonthier,  2010;  Gonthier,  2010;  Chakravarthy  and  Gonthier, 

2011;  Chakravarthy,  et  al.,  2013  A;  Chakravarthy,  2014]  Figure  3  illustrates  the  predicted 
spatial  variation  in  effective  component  pressure  Pi  ( i  =  e,  m)  and  Von-Mises  stress  t1  for 
strength  and  pressure  dominated  shocks  corresponding  to  Up  =  100  m/s  and  500  m/s,  respec¬ 
tively;  the  waves  are  propagating  from  left  to  right  in  the  figure.  Also  shown  for  comparison 
are  predicted  profiles  for  neat  HMX  corresponding  to  the  same  piston  speeds.  Spatial  vari¬ 
ations  in  the  standard  deviation  (a)  of  the  fluctuations  are  shown  to  highlight  their  local 
intensity  about  the  effective  states. 
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Figure  3:  Predictions  for  effective  component-specific  spatial  profiles  within  a  (a)  strength  dominated 
shock  ( Up  =  100  m/s)  and  (b)  pressure  dominated  shock  (Up  =  500  m/s).  The  metal  mass  fraction 
is  Xm  =  0.20.  The  effective  phase-specific  pressure  and  Von-Mises  stress  is  shown  in  each  case. 
Standard  deviations  about  the  effective  profiles  are  shown  as  thin  dashed  curves.  Black  curves 
correspond  to  neat  HMX. 


The  metal  pressure  within  the  strength  dominated  shock  is  marginally  higher  than  that  of  the 
explosive,  and  the  intensity  of  its  fluctuations  is  considerably  higher.  Moreover,  fluctuations 
in  metal  pressure  appreciably  exceed  the  effective  value  behind  the  shock  [Pm(a) / Pm  ~  2.4] 
which  is  indicative  of  high  amplitude  (relative  to  the  effective  pressure) ,  high  frequency  acous¬ 
tic  ringing  within  the  metal.  The  effective  Von-Mises  stress  for  the  metal  and  explosive  rnono- 
tonically  increase  within  the  shock  to  values  lower  than  their  corresponding  yield  strength 
which  is  indicative  of  a  small  amount  of  inelastic  deformation  at  the  particle  scale.  The 
intensities  of  Von  Mises  stress  fluctuations  are  low  relative  to  their  effective  values.  For  ap¬ 
proximately  Xm  <  0.3,  the  compaction  dynamics  is  dominated  by  explosive-explosive  contact; 
therefore,  only  small  differences  between  the  explosive  component  and  the  neat  HMX  re¬ 
sponses  are  predicted.  Though  not  shown  in  the  figure,  the  effective  dissipative  work  within 
the  metal  is  larger  than  that  of  the  explosive  due  to  its  lower  yield  strength,  and  the  in- 
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tensities  of  the  fluctuations  are  high  relative  to  their  averages,  particularly  for  the  explosive 
component  \W  d,e{&) /W  d,e  ~  2.5]  due  to  highly  localized  deformation. 

Unlike  the  strength  dominated  shock,  appreciably  higher  explosive  component  pressure  is 
predicted  behind  the  pressure  dominated  shock  than  for  neat  HMX  due  to  stronger  wave 
interactions  occurring  in  the  vicinity  of  explosive-metal  contact  surfaces.  The  high  metal 
acoustic  impedance  causes  the  amplification  of  compression  waves  within  explosive  particles 
when  reflected  from  metal  interfaces  resulting  in  a  net  increase  in  explosive  pressure.  Again, 
both  the  effective  metal  pressure  and  its  fluctuations  are  higher  than  those  of  the  explosive. 
Fluctuations  in  both  explosive  and  metal  pressure  are  appreciably  lower  than  their  effective 
values  behind  the  shock  [Pe(a)/Pe  ~  0.18  and  Pm{&)/Pm  ~  0.52]  indicating  more  spatially 
homogeneous  fields,  though  substantial  acoustic  ringing  still  occurs  within  metal  particles. 
The  effective  explosive  and  metal  Von-Mises  stresses  overshoot  their  respective  yield  surfaces 
(u>,e  =  0-37  GPa,  T0im  =  0.25  GPa)  before  relaxing  to  quasi-steady  values  behind  the  shock. 
The  Von-Mises  stresses  are  sufficiently  large  to  cause  complete  material  consolidation,  where 
compaction  of  explosive  particles  is  primarily  responsible  for  the  reduction  in  porosity  because 
the  mass  fraction  of  metal  is  low.  The  explosive  dissipative  work  is  comparable  to  that  of 
the  metal  in  this  case,  and  is  larger  than  that  of  neat  explosive  due  to  pressure  enhanced 
compaction. 

3.  Hot-spot  mass  fraction  curves  within  and  behind  compaction  shocks  were  com¬ 
putationally  analyzed.  [Chakravarthy,  et.  al.,  2010;  Chakravarthy  and  Gonthier,  2010; 
Gonthier,  2010;  Chakravarthy  and  Gonthier,  2011;  Chakravarthy,  et  ah,  2013  A]  Consider¬ 
able  attention  was  placed  in  this  effort  on  characterizing  shock  profiles  in  terms  of  compaction 
induced  dissipation  rate  because  it  is  it  believed  to  be  a  key  macroscale  hot-spot  formation 
mechanism.  Whether  the  dissipation  rate  is  locally  sufficient  to  overcome  conductive  cooling 
and  promptly  ignite  the  material  for  a  given  problem  depends,  in  part,  on  its  intensity  and 
shock  rise  time.  Therefore,  it  was  instructive  to  characterize  the  spatial  and  temporal  varia¬ 
tion  in  integrated  hot-spot  mass  within  and  behind  compaction  shocks;  the  resulting  profiles 
are  referred  to  as  hot-spot  mass  fraction  curves. 

Figure  4  illustrates  the  spatial  and  temporal  evolution  of  explosive  component  hot-spot  mass 
fraction  for  the  strength  and  pressure  dominated  shocks  discussed  above.  Color  bars  in 
these  plots  represent  the  logarithm  of  mass-fraction  (i.e. ,  log  mi)  which  better  highlights  the 
small  amounts  of  mass  that  are  locally  heated  to  elevated  temperature.  Also  shown  in  these 
plots  is  the  effective  component-specific  dissipative  heating  rate  to  indicate  the  location  and 
width  of  the  shock.  The  fraction  of  heated  mass  steadily  increases  with  distance  behind  the 
shock  until  a  quasi-steady  distribution  is  attained  near  the  end  of  the  shock  deformation 
zone.  The  structure  of  explosive  hot-spot  mass  fraction  curves  generally  consists  of  three 
parts  based  on  temperature  rise.  The  first  part  includes  mass-fractions  within  the  range 
10”1  <  me  <  10-3  that  experience  a  temperature  rise  of  AT  <  AT&  ~  100  K;  this  heated 
mass  is  almost  exclusively  due  to  volumetric  dissipation  by  plastic  work.  The  second  part 
of  the  distribution  includes  mass-fractions  within  the  range  10-3  <  me  <  10-4  that  are 
heated  through  a  temperature  rise  of  AT),  <  AT  <  A Tf  ~  200  K;  this  heated  mass  is 
due  to  both  volumetric  and  surface  dissipation,  and  it  is  stochastically  dependent  on  the 
microstructure.  The  third  part  (or  tail)  of  the  distribution  includes  mass-fractions  of  me  ~ 
10~5  that  experience  AT  >  AT)-;  this  heated  mass  is  mostly  due  to  surface  dissipation  by 
friction  work,  and  is  close  to  the  resolution  limit  of  the  computations  determined  by  the 
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Figure  4:  Predictions  for  hot-spot  mass  fraction  curves  within  a  (a)  strength  dominated  shock 
{Up  =  100  m/s)  and  (b)  pressure  dominated  shock  (Up  =  500  m/s).  The  metal  mass  fraction  is 
Xm  =  0.20.  The  black  curves  are  the  effective  dissipative  heating  rates  of  the  (a)  explosive  and  (b) 
metal  components. 


mass  of  a  single  finite-element.  The  discrete  appearance  of  the  distribution  tail  is  due  to 
numerical  discretization  of  the  temperature  field.  For  these  computations,  thermal  diffusion 
lengths  for  HMX  are  smaller  than  the  finite  element  size  (Le/ement/ Ldif fusion  ~  10).  Though 
thermal  diffusion  is  included  in  the  model,  it  has  no  discernible  effect  on  the  tail  of  the 
explosive  component  hot-spot  curve;  as  such,  explosive  hot-spots  are  effectively  adiabatic 
over  simulation  time-scales.  The  metal  exhibits  trends  that  are  qualitatively  similar  to  the 
explosive  component  with  the  exception  that  the  high  metal  thermal  diffusivity  suppresses 
the  formation  of  small,  high  frequency,  frictionally  induced  hot-spots  within  the  metal  behind 
the  shock.  Unlike  the  explosive  component,  thermal  conduction  within  the  metal  is  resolved 
on  the  computational  grid  ( Leiement/Ldiffusion  «  0.1). 

Shock  thickness  (and  particle  residence  time  within  the  shock)  significantly  decreases  for  pres¬ 
sure  dominated  shocks.  As  with  strength  dominated  shocks,  predicted  hot-spot  mass  fraction 
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curves  for  the  metal  differ  significantly  from  that  of  the  explosive  due  to  the  absence  of  a 
frictionally  induced  high-temperature  tail.  Moreover,  all  metal  mass  experiences  a  temper¬ 
ature  rise  of  at  least  A Tm  ~  80  K  and  thermal  diffusion  causes  a  slight  reduction  in  metal 
temperature  with  increasing  distance  behind  the  shock  leading  to  more  spatially  homoge¬ 
neous  fields  within  particles.  Again,  only  slight  differences  between  the  explosive  component 
heating  response  and  that  of  neat  explosive  are  predicted  because  the  fraction  of  metal  mass 
is  low.  The  temperature  rises  predicted  here  are  sufficient  to  trigger  local  adiabatic  ignition 
of  small  amounts  of  explosive  mass,  though  it  remains  unclear  whether  a  measurable  global 
ignition  event  would  subsequently  result. 

4.  The  dependence  of  shock  profiles  on  A1  mass  fraction  was  examined.  [Chakravarthy, 
et.  al.,  2010;  Chakravarthy  and  Gonthier,  2010;  Gonthier,  2010;  Chakravarthy  and  Gonthier, 
2011;  Chakravarthy,  et  al.,  2013  A;  Chakravarthy,  2014]  As  mentioned,  an  increase  in  Al 
content  increases  the  effective  shock  pressure  of  HMX-A1  composites  for  a  prescribed  piston 
speed  due  to  an  increase  in  acoustic  impedance  of  the  composite.  Though  predicted  pressure 
variations  are  qualitatively  similar  for  strength  and  pressure  dominated  shocks,  more  appre¬ 
ciable  increases  in  component  pressures  are  predicted  behind  pressure  dominated  shocks  as 
Xm  increases.  However,  contrary  to  strength  dominated  shocks,  the  effective  plastic  work  in 
both  the  explosive  and  metal  are  predicted  to  increase  with  Xm  for  pressure  dominated  shocks 
which  is  indicative  of  pressure  enhanced  compaction  followed  by  compression  of  consolidated 
material. 

Because  the  inclusion  of  metal  enhances  the  effective  explosive  dissipative  work  by  pressure 
dominated  shocks,  it  increases  explosive  hot-spot  mass.  Moreover,  an  increase  in  Xm  causes  a 
small  increase  in  the  fraction  of  explosive  mass  affected  by  plastic  work  and  slightly  suppresses 
the  high  temperature  tail  affected  by  friction  work.  The  increase  in  heating  by  plastic  work 
is  consistent  with  pressure  enhanced  compaction.  It  is  possible  that  reactive  hot-spots  which 
originate  in  the  vicinity  of  metal  particles  may  be  intensified  by  exothermic  metal  oxidation 
reactions,  thereby  facilitating  transition  to  global  ignition  by  hot-spot  growth  and  coalescence, 
provided  that  the  heat  generation  rate  is  sufficient  to  overcome  thermal  conduction  losses  and 
prevent  quenching. 

5.  The  dependence  of  shock  profiles  on  the  local  friction  coefficient  was  examined. 

[Chakravarthy,  et.  al.,  2010;  Chakravarthy  and  Gonthier,  2010;  Gonthier,  2010;  Chakravarthy 
and  Gonthier,  2011;  Chakravarthy,  et  al.,  2013  A;  Chakravarthy,  2014]  Friction  work  occurring 
in  the  vicinity  of  interparticle  contacts  is  an  important  mechanism  for  the  formation  of  small 
but  intense  hot-spots  that  result  in  local  ignition.  Though  several  studies  have  investigated 
the  role  of  frictional  heating  in  the  ignition  of  energetic  solids  based  on  Coulomb-type  models, 
disparities  exist  in  the  reported  value  of  friction  coefficient.  These  values  vary  over  the  range 
0.1  <  [i  <  1.1  for  mixtures  containing  common  explosives  such  as  HMX,  PETN,  RDX,  and 
TNT  that  are  in  sliding  contact  with  boundaries.  Specifying  an  appropriate  value  of  friction 
coefficient  is  complex,  and  depends  on  the  thermomechanical  properties  of  the  mixture  and 
boundary,  the  mixture-boundary  slip  velocity,  the  applied  normal  stress,  and  the  temperature. 
Because  of  this  complexity,  and  because  many  of  these  studies  focus  on  frictional  heating  along 
macroscale  boundaries  rather  than  local  interparticle  contacts,  it  is  desirable  to  examine  how 
the  value  of  friction  coefficient  affects  bulk  wave  propagation  and  hot-spot  formation  within 
large  particle  ensembles.  To  this  end,  predictions  were  obtained  for  several  values  of  friction 
coefficient  within  the  range  0  <  /i  <  1  for  Up  =  500  m/s  and  Xm  =  0.2.  For  simplicity,  the 
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value  of  /x  was  taken  to  be  constant  for  a  given  simulation  and  independent  of  the  type  of 
interparticle  contact  (e.g.,  HMX-HMX,  metal- metal,  or  HMX- metal). 

Predicted  Hugoniots  in  the  Up-D  and  Up-P  planes  vary  with  /x.  Uniaxial  wave  speed  D 
monotonically  increases  with  //  for  strength  dominated  shocks  as  the  rigidity  of  the  ensemble 
is  enhanced  by  interparticle  friction,  whereas  friction  only  minimally  affects  D  for  pressure 
dominated  shocks.  The  relation  D  =  D{UP)  is  approximately  linear  for  frictionless  contact 
(/x  =  0)  suggesting  that  phenomena  or  microstructural  features  that  enhance  particle  mobility 
may  give  rise  to  a  similar  result;  such  phenomena/features  might  include  particle  fracture 
with  subsequent  rearrangement,  3D  microstructures,  or  microstructures  having  higher  initial 
porosity.  The  effective  pressure  is  found  to  be  insensitive  to  variations  in  /x  for  both  strength 
and  pressure  dominated  shocks. 

Explosive  hot-spot  mass  within  and  behind  pressure  dominated  shocks  also  varies  with  /x.  The 
high  temperature  tail  of  the  explosive  hot-spot  mass  fraction  curve  is  absent  for  frictionless 
contact  (/x  =  0).  As  the  value  of  friction  coefficient  increases,  high  temperature  heating 
quickly  intensifies  until  a  value  of  /x  ~  0.1  is  reached;  further  increases  in  the  value  of  /x  over 
the  approximate  range  0.1  <  /x  <  0.6  have  little  additional  effect  on  the  high  temperature 
tail.  Values  of  /x  >  0.6,  which  result  in  frictionally  rigid  ensembles,  reduce  interparticle 
slip  causing  a  slight  reduction  in  the  intensity  of  high  temperature  mass.  This  non- monotonic 
variation  in  frictionally  induced  heating  with  /x  indicates  the  importance  of  properly  describing 
surface  dissipation  as  a  hot-spot  formation  mechanism  responsible  for  local  ignition.  Most 
mass  affected  by  plastic  work  is  largely  insensitive  to  variations  in  /x,  though  small  amounts 
of  mass  (e.g.,  A  ~  1CU3)  are  heated  to  marginally  higher  temperatures  because  tangential 
(shear)  stresses  are  more  effectively  transmitted  between  particles  as  n  increases  causing 
frictionally  enhanced  plasticity.  Qualitatively  similar  predictions  were  obtained  for  all  shock 
strengths  considered  in  this  effort. 


3  Analysis  of  Hot-Spots  in  Low  Density  Granular  Explosives 

3.1  Significance 

Various  micromechanical  mechanisms  of  shock  induced  hot-spot  formation  have  been  postulated. 
Representative  mechanisms  include  frictional  heating  due  to  sliding  contact,  shear  band  formation, 
rapid  compressive  heating  of  gas  contained  within  voids  during  void  collapse,  and  visco-plastic 
heating  of  solid  material  during  void  collapse.  Which  mechanism(s)  dominates  will  depend  on 
the  material  microstructure,  the  thermo-mechanical  properties  of  the  material  constituents,  and 
the  loading  conditions.  Moreover,  significant  fluctuations  in  stress  and  temperature  fields  behind 
shocks  due  to  material  heterogeneity  can  give  rise  to  complex  hot-spot  morphologies.  Here,  hot-spot 
morphology  refers  to  geometrical  characteristics  such  as  size  and  shape  which,  along  with  intensity, 
establishes  their  criticality  (i.e.,  the  likelihood  of  locally  experiencing  thermal  explosion).  Because 
most  analyses  typically  assume  a  simple  pore  geometry  and  a  single  formation  mechanism,  or  be¬ 
cause  they  focus  on  the  shock  response  of  a  single  hot-spot,  they  do  not  provide  detailed  statistical 
information  on  the  formation  of  hot-spots  behind  shocks.  Mesoscale  M&S  can  principally  account 
for  hot-spot  formation  and  provide  statistically  meaningful  information  on  hot-spot  intensity,  size, 
and  number  that  is  important  for  ignition  and  flame  spread  [3,  1,  11].  These  simulations  com¬ 
putationally  track  the  nonlinear  deformation  and  motion  of  individual  particles,  and  interactions 
between  particles,  enabling  the  spatial  and  temporal  evolution  of  hot-spot  fields  behind  shocks  and 
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their  combustion  implications  to  be  assessed.  This  effort  represents  a  fundamental  advancement 
in  the  characterization  of  how  material  microstructure  can  influence  fine-scale  dissipative  heating 
within  compaction  shocks  in  low  density  explosives. 

3.2  Goals  and  Methods 

The  key  goal  of  this  effort  was  to  computationally  examine  how  initial  porosity,  metal  mass  frac¬ 
tion,  and  metal  particle  size  affect  shock  induced  hot-spot  formation.  As  shown  in  the  previous 
section,  the  initial  porosity  of  granular  explosives  has  a  leading-order  effect  on  sensitivity,  but  the 
dissipative  mechanisms  responsible  for  hot-spot  formation  and  their  dependence  on  microstructure, 
composition,  and  shock  strength  remain  vague.  Importantly,  the  intensity,  size,  and  spatial  prox¬ 
imity  of  hot-spots,  and  the  degree  of  mechanical  and  thermal  confinement,  will  likely  control  the 
rate  of  hot-spot  growth  and  coalescence;  therefore,  these  features  represent  precursors  to  the  onset 
of  vigorous  combustion. 

The  computational  model  used  for  mesoscale  M&S  of  compaction  shocks  in  granular  explo¬ 
sives  was  described  in  Section  2.2.  Because  a  key  objective  of  this  effort  was  to  characterize 
hot-spots,  a  thresholding  technique  was  developed  to  identify  them  based  on  computed  temper¬ 
ature  fields  within  particles  behind  shocks.  To  this  end,  hot-spot  material  was  selected  by  per¬ 
forming  a  level-cut  through  T(x,  t)  at  a  predefined  threshold  temperature  Tt/t;  the  term  hot-spot 
refers  to  material  having  T  >  Tth ■  The  resulting  level-cut,  illustrated  in  Fig.  5,  produces  a  hot¬ 
spot  temperature  field  that  contains  N  hot-spots,  A1  fa  t),  M2(x,  t), . . .  ,ANfat),  enclosed  by  con¬ 
tours  T1(x,  t),  T2(x,  t), . . . ,  Tw(x,  t).  This  thresholding  technique  partitions  the  domain  Afat), 
into  two  disjoint  sets,  Ansfat)  and  A  Hsfat),  such  that  Ah  s  fat)  U  Ansfat)  =  A  fat)  and 
AHsfa  t)  PI  AHsfa  t)  =  0,  where  Ahs{x,  t)  =  A1(x,  t)  U  M2(x,  t)  U  . . .  U  ANfa  t ).  The  region 
Ansfa  t)  contains  all  hot-spot  material  in  M(x,  t)  for  a  given  material  realization  and  piston  speed, 
while  the  region  Ansfat)  contains  lower  temperature  material.  This  technique  filters  out  cooler 
material  allowing  for  a  quantitative  description  of  hot-spot  features.  It  is  noted  that  hot-spot  mor¬ 
phology  is  sensitive  to  simulation  resolution  and  the  choice  of  Tth-  In  this  study,  Tth  =  500  K  was 
chosen  which  is  slightly  above  the  /3-5  phase  transition  temperature  of  HMX;  as  such,  this  value 
reasonably  indicates  the  onset  of  chemical  activity,  though  exothermic  combustion  requires  T  >  Tth. 
and  is  dependent  on  hot-spot  size. 

A  complete  statistical  description  of  hot-spots  requires  that  the  hot-spot  temperature  field 
be  analyzed  for  every  realization  of  a  material  and  subsequently  averaged  over  all  realizations. 
Because  this  requirement  is  computationally  expensive,  ergodicity  was  assumed  which  stipulates 
that  averaging  over  all  realizations  of  an  ensemble  is  equivalent  to  averaging  over  the  volume  of 
a  single  realization  in  the  infinite  volume  limit  [16].  Though  computational  domains  used  in  this 
study  were  finite,  they  were  sufficiently  large  so  that  this  assumption  is  reasonable  in  establishing 
the  statistical  features  of  hot-spots.  Computational  experiments  performed  for  different  realizations 
of  a  material  having  (f>s0  ~  0.7  confirmed  this  assertion.  Detailed  expositions  of  this  work  can  be 
found  in  the  cited  references. 

3.3  Work  Accomplished 

Hot-spots  induced  by  compaction  shocks  were  computationally  examined  for  both  HMX  and  HMX- 
A1  composites.  Also,  a  preliminary  analysis  was  performed  to  explore  how  hot-spot  formation  rates 
induced  by  shocks  can  be  used  to  establish  an  ignition  criterion  for  low  density  HMX.  Such  a 
criterion  can  be  used  in  the  development  of  microstructure-dependent  ignition  and  burn  models 
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Figure  5:  Illustration  of  the  temperature  thresholding  technique  used  to  identify  hot-spots.  Here, 
Tth  =  500  K  denotes  the  temperature  threshold. 


for  assessing  transition  to  detonation,  as  addressed  in  Section  4  of  this  report.  Work  accomplished 
as  part  of  this  effort  along  with  key  observations  are  briefly  summarized  below.  Publications  and 
presentations  resulting  from  this  work  are  indicated  by  the  cited  references. 

1.  The  influence  of  particle  packing  density  on  hot-spot  formation  in  low  density 
HMX  was  computationally  examined.  [Gilbert,  2012;  Gilbert  and  Gonthier,  2012; 
Chakravarthy,  et  ah,  2013  B;  Gilbert,  et.  ah,  2013]  Emphasis  was  placed  on  examining  how  the 
material’s  initial  particle  packing  density,  characterized  by  its  effective  solid  volume  fraction 
( ps  0 ,  affects  hot-spot  statistics  for  pressure  dominated  waves  corresponding  to  piston  speeds 
within  the  range  300  <  Up  <  500  m/s.  Figure  6  shows  representative  initial  particle  ensem¬ 
bles  simulated  in  this  effort.  Ensembles  consisted  of  hexagonally  and/or  circularly  shaped, 
randomly  packed,  HMX  particles  having  diameters  of  40  /i m,  60  /jrn.  and  80  /im  in  equal 
ratios.  The  distribution  in  solid  volume  fraction  within  each  ensemble  was  nearly  Gaussian, 
and  its  standard  deviation  decreased  with  increasing  density. 

Analysis  of  temperature  field  predictions  indicate  that  hot-spot  intensity  is  only  marginally 
affected  by  changes  in  impact  speed,  initial  porosity,  and  particle  shape  for  0.57  <  <^s  0  <  0.84, 
with  higher  temperature  scatter  existing  for  highly  porous  material  due  to  the  wider  distri¬ 
bution  of  initial  pore  sizes.  Unlike  intensity,  hot-spot  morphology  and  spatial  proximity  vary 
more  appreciably.  Figure  7  summarizes  the  predicted  variation  in  hot-spot  number  density, 
volume  fraction,  and  specific  surface  area  with  Up  for  each  ensemble.  These  quantities,  ex¬ 
pressed  per  unit  total  volume  of  material  behind  the  shock,  are  often  explicitly  or  implicitly 
used  to  describe  ignition  phenomena  in  macroscale  burn  models.  An  exponential  growth  in 
each  of  these  quantities  is  predicted  over  the  range  of  piston  speeds  considered.  Importantly, 
these  quantities  exhibit  a  nonlinear  dependence  on  <^s0:  packing  densities  within  the  range 
0.68  <  4>s  o  <  0.84  result  in  substantially  fewer  hot-spots  at  a  given  impact  speed,  whereas  no 
appreciable  changes  are  predicted  for  approximately  0.68  <  4>s  0  <  0.58.  This  latter  predic¬ 
tion  is  consistent  with  shock  sensitivity  data  for  neat  granular  explosives  which  indicate  little 
variation  in  sensitivity  with  porosity  for  highly  porous  materials.  The  slopes  of  curves  cor¬ 
responding  to  different  ensembles  are  similar,  with  the  exception  of  hot-spot  number  density 
which  has  a  lower  rate  of  increase  with  Up  for  highly  porous  material  due  to  the  clustering 
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Figure  6:  Representative  initial  particle  ensembles  used  with  mesoscale  M&S  to  assess  the  effect  of 
packing  density  on  hot-spot  formation  in  low  density  HMX. 
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Figure  7:  Predicted  variation  in  (a)  hot-spot  number  density,  (b)  hot-spot  volume  fraction,  and 
(c)  hot-spot  specific  surface  area  for  each  ensemble  with  piston  speed.  The  dashed  lines  represent 
best-fit  exponential  functions. 


and  coalescence  of  hot-spots  at  higher  impact  speeds;  i.e. ,  hot-spots  coalesce  as  more  material 
is  heated  to  elevated  temperature,  and  their  number  density  decreases.  This  observation  is 
further  addressed  in  Item  2  of  this  section. 

The  effective  volumetric  strain  ey  =  1  —  p0/p  is  influenced  by  both  Up  and  (j)s0 ,  where 
p  is  the  effective  material  density.  For  pressure  dominated  shocks  having  </>  ~  1  with  a 
stiff  solid  equation  of  state  so  that  solid  density  changes  are  small,  the  effective  volumetric 
strain  may  be  approximated  by  ey  ~  l  —  cf)s  0/cj)s  ~  1  —  (f>s  0  which  is  the  initial  porosity. 
Predicted  deformation  and  temperature  fields  appreciably  differ  depending  on  (f)s0  with  high 
volumetric  deformation  resulting  in  larger,  more  densely  packed  hot-spots  having  a  higher 
surface  area/volume  ratio.  Consequently,  the  gas  phase  production  rate  responsible  for  the 
onset  of  vigorous  combustion  behind  sustained  shocks  is  more  likely  to  be  controlled  by 
hot-spot  morphology,  proximity,  and  number  rather  than  temperature-dependent  chemistry. 
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Figure  8:  An  illustration  of  hot-spot  axial  spacing  used  to  compute  the  formation  time  between 
successive  hot-spots  based  on  their  centroid  locations. 


Though  speculative  in  the  absence  of  coupled  multiphase  transport  and  chemistry,  a  similar 
description  has  been  postulated  to  explain  shock  sensitivity  data  for  preheated  TATB-based 
explosives.  These  results  suggest  that  the  observed  increase  in  sensitivity  with  initial  porosity 
for  sustained  loading  is  likely  due  to  an  increase  in  hot-spot  size  and  number  rather  than 
intensity. 

2.  The  dynamics  of  hot-spots  in  low  density  HMX  was  characterized  in  terms  of  in¬ 
creasing  shock  strength.  [Chakravarthy,  2014;  Gonthier  and  Chakravarthy,  2014;  Rao  and 
Gonthier,  2015  A]  This  analysis  focused  on  describing  how  hot-spots  evolve  with  compaction 
shock  strength;  therefore,  it  was  convenient  to  analyze  hot-spots  in  a  wave-attached  frame, 
as  conceptually  shown  in  Fig.  8.  In  this  frame,  unstressed  porous  material  enters  the  wave  at 
speed  D  and  stressed  compacted  material  leaves  the  wave  at  speed  w  =  D  —  Up,  where  Up  is 
the  piston  speed  and  D  is  the  corresponding  steady  wave  speed.  The  spatial  shock  structure 
has  a  finite  thickness  5  due  to  energy  dissipation  and  dispersion  within  the  wave.  Though 
the  compacted  state  possesses  fluctuations  in  thermomechanical  fields,  its  effective  pressure 
or  temperature  can  be  used  to  characterize  shock  strength.  Hot-spots  are  formed  within  the 
shock  due  to  inelastic  pore  collapse,  emerging  from  the  shock  at  a  frequency  that  depends  on 
both  the  mesostructure  and  shock  strength. 

Individual  hot-spots  are  identified  by  the  location  of  their  area  centroids  enabling  the  for¬ 
mation  rate  to  be  characterized  based  on  hot-spot  spacing.  Though  other  metrics  may  be 
used  to  characterize  spacing  (e.g.,  surface-to-surface  nearest  neighbor  distributions,  point-to- 
surface  distributions,  etc.),  numerical  experiments  indicate  that  the  formation  rate  is  largely 
insensitive  to  the  choice  of  metric  for  shocks  dominated  by  hot-spot  nucleation,  as  discussed 
below.  For  a  given  simulation,  the  formation  time,  defined  as  the  time  between  successive 
formation  of  hot-spots  within  the  shock,  is  estimated  by  r/  =  5hs/(D  ~  Up),  where  Shs  is 
the  axial  distance  between  their  centroids  in  the  compacted  state  identified  in  Fig.  8.  The 
transverse  distance  between  hot-spots  is  not  used  to  compute  the  distribution  of  formation 
times,  but  is  implicitly  accounted  for  by  expressing  the  effective  formation  rate  as  an  intensive 
quantity  using  the  transverse  wave  dimension  d  in  the  analysis  that  follows.1  Development 

1The  wave  dimension  d  is  given  by  the  transverse  width  of  the  computational  domain. 
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Figure  9:  Predicted  variation  in  hot-spot  formation  time  distribution  with  Up  for  steady  waves  in 
porous  HMX:  (a)  4>0  =  0.68;  (b)  4>0  =  0.84. 


of  a  technique  that  explicitly  accounts  for  the  average  transverse  spacing  between  hot-spots 
is  a  current  topic  of  our  ongoing  research.  This  definition  of  tj,  which  establishes  the  effec¬ 
tive  formation  rate,  is  influenced  by  both  hot-spot  spacing  and  wave  speed.  It  is  plausible 
that  mesostructures  having  comparable  hot-spot  spacing  but  different  wave  speeds  will  have 
different  ignition  times  because  their  formation  rates  will  differ. 

Figure  9  illustrates  predicted  distributions  in  hot-spot  formation  time  for  steady  shocks  in 
porous  HMX.  The  distribution  is  dependent  on  both  shock  strength  and  mesostructure,  char¬ 
acterized  by  (j) q  =  0.68  and  0.84.  Because  simulations  of  weak  impact  require  long  com¬ 
putational  domains  to  obtain  statistically  significant  numbers  of  sparsely  distributed  hot¬ 
spots,  this  analysis  focuses  on  shocks  corresponding  to  Up  >  400  m/s  which  have  sufficient 
strength  to  eliminate  porosity.  Formation  time  distributions  are  approximately  exponential 
indicating  that  hot-spot  formation  represents  a  time-invariant  Poisson  process  for  mesostruc¬ 
tures  considered  in  this  effort.  The  cumulative  exponential  distribution  function  is  given 
by  F(rf)  =  1  —  exp  AryJ,  where  A  is  the  effective  formation  rate/frequency  of  hot-spots. 
Because  the  mechanics  of  hot-spot  formation  depend  on  the  interparticle  contact  geometry, 
hot-spot  formation  times  may  be  strongly  influenced  by  the  particle  size  distribution. 

Predicted  variations  in  A  with  shock  strength  and  mesostructure,  summarized  in  Fig.  10  in 
terms  of  Up  for  </>0  =  0.68,  0.77,  and  0.84,  can  provide  insight  into  the  dynamics  of  hot¬ 
spot  formation  within  shocks  that  may  have  ignition  implications.  The  dynamics  involves  a 
competition  between  the  nucleation,  growth,  and  agglomeration  of  hot-spots  as  shock  strength 
increases.  Here,  nucleation  refers  to  the  formation  or  seeding  of  new  hot-spots  in  the  vicinity 
of  interparticle  contact  surfaces,  whereas  growth  and  agglomeration  refer  to  an  increase  in 
hot-spot  size  and  the  coalescence  of  neighboring  hot-spots  due  to  uncontained  plastic  flow, 
respectively.2  For  each  material,  an  exponential  increase  in  A  is  initially  predicted  with 
increasing  Up  due  to  the  influence  of  hot-spot  nucleation,  with  higher  porosity  resulting  in 
a  higher  formation  rate  for  fixed  Up.  As  Up  further  increases,  neighboring  hot-spots  begin 

J  Because  reactive  hot-spots  are  identified  by  their  area  centroids  in  this  analysis  which  is  independent  of  their 
size,  agglomeration  can  cause  in  a  reduction  in  A  when  it  dominates  nucleation. 
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Figure  10:  Predicted  variation  in  the  effective  hot-spot  formation  rate  with  (a)  piston  speed  and 
(b)  shock  pressure  for  low  density  HMX. 


to  coalesce/agglomerate  resulting  in  slower  increases  in  A;  the  value  of  Up  associated  with 
the  onset  of  this  “strong  shock”  regime  increases  with  4>0.  In  this  regime,  the  value  of  A 
may  decrease  with  a  further  increase  in  wave  strength  as  agglomeration  begins  to  dominate 
nucleation  and  hot-spots  become  ubiquitous  within  the  material;  this  scenario  is  illustrated 
by  the  case  corresponding  to  Up  =  900  m/s  and  (j>0  =  0.84.  It  is  plausible  that  waves  within 
the  strong  shock  regime  will  result  in  prompt  ignition  and  vigorous  burn  of  material  due  to 
the  close  spatial  proximity  of  large  and  intense  hot-spots. 

3.  The  effective  hot-spot  formation  rate  in  low  density  HMX  can  be  expressed  by 
a  power-law  relation  in  terms  of  shock  pressure.  [Chakravarthy,  2014;  Gonthier  and 
Chakravarthy,  2014;  Rao  and  Gonthier,  2015  A]  It  is  desirable  to  express  the  formation  rate 
in  intensive  form  based  on  the  wave  frontal  area.  For  the  2-D  simulations  performed  in 
this  study,  the  specific  formation  rate  was  defined  by  A  =  X/d,  where  d  is  the  transverse 
wave  width.  Therefore,  A  represents  the  hot-spot  number  flux  emerging  from  the  wave.  The 
formation  rate  was  then  related  to  shock  pressure  by  the  power-law  expression  A  =  aP  . 
The  data  and  best-fit  curves  for  each  </>0  are  shown  in  Fig.  10  for  the  nucleation  phase.  Both 
the  prefactor  a  and  exponent  n  vary  with  initial  packing  density.  As  </0  increases  for  fixed 
P,  the  effective  bulk  modulus  of  the  material  increases  resulting  in  volumetric  stiffening  and 
faster  waves,  and  reactive  hot-spot  spacing  increases  due  to  enhanced  stress  transmission 
between  particles;  consequently,  A  decreases.  However,  because  n  =  Slog  X/d  log  P,  the 
reactive  hot-spot  formation  rate  increases  more  rapidly  with  pressure  for  dense  material, 
though  it  remains  relatively  low  throughout  the  nucleation  regime.  This  more  rapid  increase 
in  formation  rate  with  pressure  for  dense  material  is  consistent  with  abrupt  ignition  once  a 
critical  threshold  is  obtained.  These  predictions  indicate  that  the  effective  hot-spot  formation 
rate  is  sensitive  to  material  mesostructure  and,  therefore,  may  be  a  plausible  measure  for 
establishing  a  mesostructure-dependent  ignition  theory  for  macroscale  modeling  of  low  density 
granular  explosives.  This  topic  is  further  addressed  in  Section  4  of  this  report. 

4.  A  “Pnr-type”  relation  for  the  ignition  time  behind  a  shock  was  established  based 
on  a  cumulative  ignition  function  that  depends  on  mesostructure  through  the  ef¬ 
fective  hot-spot  formation  rate.  [Gonthier  and  Chakravarthy,  2014]  The  ignition  relation 
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assumes  that  1)  ignition  occurs  during  hot-spot  nucleation  so  that  the  ignition  time  r*  >  rr, 
where  Tr  is  the  shock  rise  time,  and  2)  r*  ~  A-1.  The  first  assumption  restricts  the  analysis 
to  shocks  of  intermediate  strength  that  have  an  induction  period  prior  to  ignition  due  to 
growth  and  coalescence  of  initially  small  hot-spots  by  multiphase  transport  mechanisms.  The 
second  assumption  is  plausible  because  a  large  value  of  A  corresponds  to  the  rapid  formation 
of  closely  packed  hot-spots  that  is  likely  to  enhance  the  onset  of  ignition.  Implicit  in  these 
assumptions  is  that  waves  resulting  in  significant  hot-spot  agglomeration  due  to  widespread 
plasticity  induce  prompt  ignition  because  t*  <  tt. 

The  cumulative  hazard  (ignition)  function  corresponding  to  an  exponential  distribution  of 
hot-spot  formation  times  for  a  shock  characterized  by  A  is  given  by  Z(r)  =  /Qr  A  dr'  =  At. 
Here,  Z(r)  represents  the  specific  accumulation  of  hot-spots  behind  the  shock  over  time.  As 
a  first  approximation,  it  is  assumed  that  ignition  occurs  behind  the  shock  when  Z  =  Z*. 
The  ignition  time  is  given  by  r*  =  Z*A~X  which  can  be  expressed  in  terms  of  P  and  4>0. 
The  relation  for  ignition  time  is  then  given  by  P"t*  =  The  value  of  Z*  is  estimated 
based  on  DDT  data  for  porous  HMX  obtained  by  McAfee,  et  al.  [10],  which  indicate  that 
t*  ~  400  /is  for  P  =  0.04  GPa,  0O  =  0.66,  and  d  =  125  /mi.  Using  these  values  with 
the  power-law  parameter  values  for  (j)0  =  0.68,  the  estimate  Z*  =  7.27  nun-1  is  obtained. 
Assuming  this  value  to  be  a  property  of  HMX,  it  is  possible  to  estimate  the  variation  in 
ignition  time  as  a  function  of  shock  strength  for  different  values  of  (f>Q.  This  assumption  is 
particularly  reasonable  for  shocks  considered  in  this  study  which  have  sufficient  strength  to 
eliminate  porosity  because  the  effect  of  residual  porosity  on  flame  spread  leading  to  global 
ignition  is  inconsequential;  therefore,  the  ignition  time  will  be  controlled  to  leading-order  by 
the  spatial  proximity  of  reactive  hot-spots  in  the  completely  compacted  state  and  the  rate  at 
which  they  are  formed. 

Figure  11  illustrates  the  predicted  variation  in  ignition  time  behind  steady  waves  in  porous 
HMX.  Highlighted  in  red  on  the  (f>0  =  0.68  curve  is  the  single  datum  used  to  establish  the  value 
of  Z*.  Also  highlighted  in  the  figure  is  the  approximate  range  of  measured  shock  rise  times 
reported  by  Sheffield,  et  al.  [12],  for  porous  HMX  for  shock  strengths  considered  in  this  study 
(0.07  <  Tr  <  0.7  /is).  Measurements  indicate  that  shock  rise  time  is  largely  insensitive  to  the 
value  of  <p0.  The  predicted  ignition  time  rapidly  decreases  with  increasing  wave  strength  for 
fixed  <f>0.  There  exists  a  critical  shock  pressure  that  results  in  r*  ~  rr  indicating  the  onset  of 
prompt  ignition  within  the  wave;  interestingly,  this  critical  pressure  approximately  equals  that 
required  for  the  onset  of  hot-spot  agglomeration.  Using  the  terminology  of  Ref.  [12],  Region 
I  identified  in  the  figure,  which  corresponds  to  t*  >  rr,  represents  an  ignition  controlled 
process  because  of  the  long  induction  times,  whereas  Region  II,  which  corresponds  to  r*  <  rr, 
represents  a  growth  controlled  process  because  the  waves  are  sufficiency  strong  to  induce 
prompt  ignition.  Measurements  indicate  that  the  reaction  zone  is  strongly  coupled  to  the 
shock  front  for  growth  controlled  processes  which  is  consistent  with  the  estimates  given  here. 
Also  shown  in  Fig.  11(b)  are  measured  detonation  initiation  times  reported  by  McAfee,  et 
al.,  [10]  and  Dick  [7]  for  <fiQ  ~  0.68.  The  relative  difference  between  the  ignition  and  initiation 
time  increases  with  shock  pressure  indicating  that  reaction  growth  following  ignition  plays  an 
increasingly  important  role  in  establishing  the  initiation  time  scale. 

The  relative  shock  sensitivity  of  explosives  is  frequently  reported  in  the  literature  in  terms  of 
the  axial  run  distance  to  detonation  as  a  function  of  input  shock  pressure.  Assuming  that  run 
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Figure  11:  (a)  Predicted  variation  in  the  global  ignition  time  with  wave  strength  for  porous  HMX 
having  different  values  of  (f>0.  (b)  Comparison  of  the  predicted  ignition  curve  with  initiation  time 
data  for  </0  =  0.68. 


distance  is  a  monotonic  function  of  ignition  time,  then  it  is  possible  to  qualitatively  compare 
measurements  of  run  distance  to  ignition  time  predictions.  In  particular,  measurements  on 
low  density  HMX  indicate  that  the  run  distance  to  detonation  increases  with  (j)0  for  fixed  shock 
pressure,  and  that  it  is  slightly  more  sensitive  to  variations  in  pressure  for  dense  material  [8, 
12].  Predicted  ignition  time  curves  reflect  both  features,  where  the  increase  in  pressure 
sensitivity  is  due  to  the  increase  in  the  value  of  power-law  exponent  n  with  <f>0.  These  results 
were  used  to  establish  an  ignition  and  burn  model  for  describing  the  shock  initiation  of  low 
density  HMX,  as  discussed  in  Section  4  of  this  report. 

5.  The  influence  of  particle  packing  density  on  hot-spot  formation  in  low  density 
HMX-A1  composites  was  computationally  examined.  [Chakravarthy,  2014;  Chakravarthy 
and  Gonthier,  2016]  Particle  ensembles  were  considered  in  this  effort  that  isolated  the  effects 
of  initial  particle  packing  density,  characterized  by  the  solid  volume  fraction  (f>s ,  and  aluminum 
mass  fraction  Am  on  shock  sensitivity.  Representative  ensembles  used  with  mesoscale  M&S 
are  shown  in  Fig.  12,  where  yellow  particles  are  HMX  and  red  particles  are  Al.  The  ensembles 
consisted  of  hexagonal  and/or  circular  particles  having  an  average  diameter  of  60  /mi  with 
a  narrow  distribution.  The  effective  packing  density  was  varied  over  the  range  0.678  <  cj)s 
<  0.835.  An  effective  metal  mass  fraction  of  Am  =  0.2  was  used  for  all  ensembles.  Impact 
speeds  considered  in  this  effort  varied  over  the  range  400  <  Up  <  800  m/s  which  resulted  in 
pressure  dominated  shocks. 

Figure  13  shows  predicted  Hugoniots.  Here  Up  and  D  are  the  piston  and  shock  speed,  respec¬ 
tively;  P  =  {Pe(t>e  +  Pm&m) / (4>e  +  <^m)  is  the  effective  equilibrium  mixture  pressure  behind 
the  shock,  where  Pe  and  Prn  are  the  explosive  and  metal  pressures,  and  (j)e  and  </>m  are  the  ex¬ 
plosive  and  metal  volume  fractions;  and  Wpe  is  the  effective  dissipative  work  in  the  explosive 
component.  Also  shown  are  experimentally  measured  Hugoniots  for  neat  HMX  [12].  Predic¬ 
tions  indicate  that  the  inclusion  of  aluminum  increases  wave  speed,  pressure,  and  explosive 
dissipation  for  all  cases  compared  to  equivalent  packing  densities  of  neat  HMX.  However, 
little  variation  in  explosive  dissipation  with  shock  pressure  is  predicted  due  to  the  increase 
in  pressure  with  piston  speed  for  HMX-A1  mixtures.  Shock  rise  times  were  found  to  decrease 
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Figure  12:  Initial  particle  ensembles  with  varying  packing  densities  4>s  and  composition  Xm.  Here, 
red  particles  represent  Al,  and  yellow  particles  represent  HMX. 


with  an  increase  in  packing  density  and/or  piston  speed  but  were  insensitive  to  variations  in 
metal  mass  fraction. 

An  increase  in  packing  density  was  found  to  increase  the  average  HMX  hot-spot  size  for  fixed 
piston  speed  with  only  marginal  variations  in  hot-spot  intensity.  The  number  density  of  hot¬ 
spots  based  on  the  total  volume/area  increases  with  piston  speed  for  all  packing  densities 
prior  to  the  onset  of  hot-spot  agglomeration,  whereas  the  volume  fraction  of  hot-spots  rnono- 
tonically  increases.  Relative  to  neat  HMX,  the  overall  number  density  of  hot-spots  is  lower 
due  to  spatial  dilution. 

When  based  on  HMX  volume  and  not  total  volume,  qualitatively  different  behavior  is  pre¬ 
dicted.  Hugoniots  were  used  to  correlate  hot-spot  number  density  and  volume  fraction  (based 
on  HMX  volume)  with  shock  pressure,  as  shown  in  Fig.  14.  Aluminized  HMX  is  predicted  to 
slightly  suppress  hot-spot  formation  relative  to  neat  HMX  for  low  pressure  shocks  associated 
with  the  nucleation  and  growth  phases  of  hot-spot  formation;  thus,  both  the  explosive  hot¬ 
spot  number  density  and  volume  fraction  marginally  decrease  with  metal  addition.  For  low 
pressure  shocks,  where  surface  dissipation  due  to  friction  and/or  localized  plasticity  are  pri¬ 
mary  heating  mechanisms,  metal-explosive  contact  surfaces  act  as  energy  sinks  due  to  the  high 
thermal  conductivity  of  aluminum  which  suppresses  hot-spot  formation.  As  shock  pressure 
increases  into  the  agglomeration  phase  of  hot-spot  formation,  uncontained  plasticity  within 
explosive  particles  becomes  a  primary  heating  mechanism  which  is  enhanced  by  intense  shock 
reflections  from  explosive- aluminum  interfaces.  Moreover,  this  effect  is  predicted  to  increase 
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Figure  13:  Predicted  variations  in  effective  shock  end  states  as  function  of  4>s  and  Xm:  ( a)P-Up ; 
(b)  D-Up-  (c)  Wpe-Up ■  and  (d)  Wpe-P. 


with  particle  packing  density.  These  predictions  suggest  that  the  inclusion  of  aluminum  in 
low  density  composites  may  delay  HMX  ignition  for  low  pressure  shocks. 

6.  The  influence  of  particle  size  on  hot-spot  formation  in  low  density  HMX-A1 
composites  was  computationally  examined.  [Chakravarthy,  2014;  Chakravarthy  and 
Gonthier,  2016]  Mesoscale  M&S  was  performed  to  characterize  how  disparities  between  the 
average  HMX  and  aluminum  particle  size  affect  explosive  hot-spot  formation.  To  this  end, 
aluminum/HMX  particle  size  ratios  of  1,  2,  3,  and  6  were  considered  for  aluminum  mass 
fractions  of  Xm  =  0.20  and  0.56,  and  a  solid  volume  fraction  of  (j>s  =  0.84.  All  particle  size 
distributions  were  narrow,  and  the  average  HMX  particle  size  was  approximately  60  /im  for 
all  cases.  These  particle  size  ratios  were  largely  chosen  based  on  computational  constraints. 
Representative  particle  ensembles  are  shown  in  Fig.  15. 

Predictions  indicate  that  aluminum  particle  size  has  a  largely  inconsequential  effect  on  shock 
Hugoniots  for  material  having  Xm  =  0.20,  as  illustrated  by  the  variation  in  explosive  dissipa¬ 
tive  work  with  shock  pressure  in  Fig.  16(a).  Compaction  shock  Hugoniots  are  more  strongly 
influenced  by  aluminum  particle  size  for  materials  having  higher  aluminum  mass,  as  illustrated 
in  Fig.  16(b)  for  Xm  =  0.56,  where  a  monotonic  increase  in  shock  pressure  and  dissipative  work 
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Figure  14:  Predicted  variations  in  the  effective  plastic  work  of  the  explosive  component  as  a  function 
of  (j)s,  and  \m. 
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Figure  15:  Meso-structures  with  varying  particle  sizes,  particle  distributions  and  composition.  Here, 
red  particles  represent  metal,  and  yellow  particles  represent  explosive. 


are  predicted  with  an  increase  in  A1  particle  size.  Because  this  size  effect  is  also  predicted  for 
identical  particle  ensembles  of  neat  HMX,  it  is  not  dominated  by  the  material  make-up  of  the 
contacting  surfaces.  Additionally,  smaller  particles  that  are  surrounded  by  larger  particles 
can  become  pinched  during  shock  passage  resulting  in  higher  local  pressure  and  deformation; 
this  pinching  effect  is  enhanced  in  regions  where  small  HMX  particles  are  locally  surrounded 
by  large  A1  particles  due  to  the  high  bulk  modulus  of  aluminum. 
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Figure  16:  Predicted  variations  in  explosive  dissipative  work  with  particle  size  ratio  and  shock 
pressure  for  a  granular  HMX-A1  composite  having  </>s  =  0.84:  (a)Am  =  0.20  and  (b)Am  =  0.56. 


Figure  17:  Predicted  variations  in  hot-spot  volumetric  properties  with  particle  size  ratio  and  shock 
pressure  for  a  granular  HMX-A1  composite  having  cj)s  =  0.84  and  Am  =  0.20:  (a)  number  density 
and  (b)  volume  fraction. 


Explosive  hot-spot  number  density  and  volume  fraction  are  predicted  to  have  a  strong  non¬ 
linear  dependence  on  A1  particle  size  for  low  pressure  shocks  (approximately  P  <  2.5  GPa  for 
4>s  =  0.84),  as  illustrated  in  Fig.  17.  An  appreciable  decrease  in  hot-spot  number  density  and 
volume  fraction  is  predicted  for  material  containing  large  A1  particles  compared  to  material 
having  equal  particle  sizes.  Though  speculative,  it  is  plausible  that  frictional  dissipation 
and/or  localized  plasticity  which  dominate  the  nucleation  and  early  growth  phase  of  hot-spot 
formation  may  be  sensitive  to  minor  changes  in  granular  bed  morphology.  Marginal  differences 
in  hot-spot  volume  fraction,  and  larger  differences  in  number  density,  are  predicted  for  high 
pressure  shocks.  This  prediction  is  due  to  uncontained  plasticity  and  hot-spot  agglomeration 
that  occurs  in  the  vicinity  of  large  A1  particles  resulting  in  the  formation  of  fewer  but  larger 
hot-spots.  Similar  trends  are  predicted  for  materials  having  Am  =  0.56.  Additional  analysis  is 
required  to  examine  why  disparate  particle  sizes  can  result  in  a  significant  decrease  in  hot-spot 
volumetric  quantities  at  low  shock  pressures;  this  effect  may  have  significant  implications  on 
the  shock  ignition  of  HMX-A1  composites. 
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4  Ignition  and  Burn  Modeling  for  Low  Density  HMX 

4.1  Significance 

Continuum  theories  are  routinely  used  to  examine  shock  initiation  of  solid  explosives  that  occurs 
over  length  scales  that  are  appreciably  larger  than  the  particle  scale.  These  theories  describe  the 
effective/bulk  response  of  the  material  and  do  not  resolve  the  formation,  growth,  and  interaction  of 
reactive  hot-spots  within  the  microstructure.  The  relative  importance  of  dissipation  mechanisms, 
their  dependence  on  microstructure  and  loading  conditions,  and  their  influence  on  macro-scale  be¬ 
havior  remain  fertile  areas  of  research.  Both  mesoscale  and  macroscale  M&S  have  been  performed 
to  address  these  and  other  issues,  but  these  studies  have  largely  focused  on  sustained  single  shocks 
to  establish  a  foundation  for  analyzing  more  complex  loading  scenarios  that  may  result  in  detona¬ 
tion.  Successful  application  of  commonly  used  hot-spot  motivated  burn  models,  such  as  Ignition 
and  Growth  [15],  hinges  on  proper  tuning  of  large  parameter  sets  to  initiation  data.  Though  many 
parameters  have  a  physical  interpretation  based  on  hot-spots,  it  can  be  difficult  to  establish  their 
values  in  the  absence  of  hot-spot  and  initiation  data  for  specific  microstructures  under  consider¬ 
ation.  This  effort  represents  a  preliminary  step  toward  the  development  of  a  simple  technique 
to  incorporate  nricrostructure-dependent  hot-spot  features  into  ignition  and  burn  models  for  low 
density  explosives  without  the  complexity  of  having  to  resolve  particle  scale  ignition  and  burn. 


4.2  Goals  and  Methods 

The  key  goal  of  this  effort  was  to  utilize  results  from  mesoscale  M&S  discussed  in  Section  3  to 
formulate  a  micro  structure- dependent  ignition  and  burn  model  for  low  density  HMX.  The  model 
was  used  with  a  macroscale  multiphase  flow  theory  and  a  high-resolution  computational  technique 
to  examine  shock  initiation.  The  model  is  conceptually  similar  to  Ignition  and  Growth  but  accounts 
for  ignition  in  terms  of  parameters  that  are  explicitly  determined  from  hot-spot  fields. 


The  multiphase  flow  theory  used  in  this  effort  is  a  two-phase  limiting  form  of  the  more  general 
multiphase  model  formulated  by  Crochet  [5].  The  vector  form  of  the  model  equations  in  Cartesian 
coordinates  are  expressed  as 


<9q  <9f(q) 

dt  dx 


g(q)^T  +  s(q)- 


(4) 


Here,  q  is  the  vector  of  conserved  variables,  f(q)  is  the  flux  vector,  g(q)  is  the  vector  associated  with 
nonconservative  nozzling  sources,  and  s(q)  is  the  phase  interaction  source  vector.  These  equations 
describe  the  evolution  of  mass,  momentum,  total  energy,  and  solid  volume  fraction  of  the  granular 
explosive.  Complete  details  of  the  model,  including  constitutive  relations  for  granular  HMX,  are 
given  in  Refs.  [5,  6].  Development  of  the  multiphase  constitutive  theory  represents  an  additional 
contribution  of  this  effort  as  per  Item  3  below.  All  mesoscale  M&S  performed  in  this  effort  was 
based  on  the  model  and  numerical  technique  discussed  in  Section  2.  Both  inert  and  reactive  shocks 
in  low  density  HMX  were  examined.  Comparisons  between  mesoscale  predictions  of  effective  shock 
profiles  and  those  given  by  the  multiphase  theory  for  inert  shocks  were  examined  to  assess  their 
agreement,  particularly  in  terms  of  dissipation;  this  assessment  was  important  because  it  enables 
particle  scale  behavior  to  be  interpreted  within  the  context  of  multiphase  theory  predictions. 


4.3  Work  Accomplished 

Work  accomplished  as  part  of  this  effort  along  with  key  observations  are  briefly  summarized  below. 
Publications  and  presentations  resulting  from  this  work  are  indicated  by  the  cited  references. 
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1.  An  analysis  of  compaction  wave  dissipation  in  low  density  HMX  was  performed. 

[Rao  and  Gonthier;  2014  A;  Rao,  et  al.,  2013;  Rao  and  Gonthier,  2016]  Mesoscale  M&S 
was  performed  to  computationally  examine  how  initial  porosity  influences  dissipation  rates 
within  inert  compaction  shocks  and  to  compare  predictions  to  a  continuum  compaction  theory. 
Effective  shock  profiles  were  obtained  by  averaging  meso-scale  fields  over  space  and  time.  The 
continuum  theory  directly  predicts  the  variation  in  effective  thermomechanical  fields  during 
compaction  due  to  an  imbalance  between  an  effective  solid  pressure  and  a  configurational 
stress. 

A  key  focus  of  this  effort  was  to  examine  dissipation  within  compaction  shocks.  The  mass- 
specific  deformation  power  prescribed  by  the  mesoscale  model  [Section  2]  is  expressed  by 

V  =  -  (r  :  (d  —  dp))  +  — tr(d  -  dp)  +  -  (cr  :  dp) .  (5) 

P  P  P 

Terms  on  the  right  side  of  this  equation  represent  shear,  compression,  and  plastic  work  rates, 
respectively.  Within  the  context  of  this  model,  shear  work  results  in  a  non-thermal  change 
of  stored  strain  energy,  whereas  compression  and  plastic  work  result  in  deformation  induced 
heating.  Only  plastic  work  is  dissipative  in  that  it  causes  an  irreversible  production  of  entropy, 
with  the  local  dissipation  power  given  by 

Trj  =  -  (er  :  dp)  >  0,  (6) 

P 

where  i]  is  the  mass-specific  entropy.  Friction  work  occurring  at  inter-particle  boundaries  is 
also  dissipative  when  slip  occurs;  the  frictional  dissipation  power  is  locally  given  by  Tiy  = 
tc  •  vr,  where  vr  is  the  relative  slip  velocity.  Because  frictional  dissipation  is  localized  at 
particle  contact  surfaces,  it  does  not  appreciably  contribute  to  the  effective  (or  bulk)  mass- 
specific  dissipated  energy.  Therefore,  dissipated  energy  by  plastic  work  dominates  the  effective 
response,  though  frictionally  induced  tangential  surface  tractions  can  indirectly  contribute  to 
dissipation  by  enhancing  plastic  deformation  within  particles.  Dissipative  work  is  obtained 
by  integrating  Eq.  (6)  in  time. 

The  mass-specific  deformation  power  prescribed  by  the  continuum  compaction  theory  is  ex¬ 
pressed  by 


where 


V 


B  +  ep  + 

JL^L  +  +  (ps  -  P)  d(j)  |  /3  d(j) 

ps(j)  dt  v  /  ps 2  dt  ps4>  dt  ps(j)  dt. 


dt  p 

d±  =  {  /W'i  if  /  (0)  >  ^ 

dt  I  0  otherwise. 


(7) 

(8) 
(9) 


Here,  d  (•)  /dt  =  d  (•)  /dt  +  v  •  V  (•)  is  the  Lagrangian  derivative.  Terms  on  the  left  of  the 
second  equality  in  Eq.  (7)  are  respectively  defined  by  the  corresponding  terms  on  the  right. 
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In  Eq.  (7),  B  implicitly  accounts  for  the  rate  of  stored  compaction  energy  at  the  particle 
scale  due  to  shear,  ep  accounts  for  the  compression  work  rate,  and  ec  =  e<f,  +  e~  accounts 
for  the  compaction  work  rate.  The  term  ec  consists  of  both  a  rate-dependent  component  e^, 
which  vanishes  in  the  slow  compaction  limit  (i.e.,  Ps  — >  /3),  and  an  inelastic  component  e~. 
These  contributions  to  the  deformation  power  are  respectively  analogous  to  the  integrated 
contributions  of  those  given  in  Eq.  (5)  for  the  mesoscale  description.  The  Second  Law  of 
Thermodynamics  gives  the  strong  form  of  the  dissipation  inequality 

TSV  =  e<f,  +  e~  >  0,  (10) 

where  p  =  r/s  is  the  mass-specific  granular  solid  entropy.  The  first  term  on  the  right  side 
of  this  equation  is  non-negative,  as  is  the  second  term  because  /3  >  0  and  dcp/dt  >  0  by 
construction.  This  thermodynamic  description  is  compatible  with  a  Helmholtz  free  energy  of 
the  form  ip(ps,  Ts,  (j>,  <j>)  =  il>s(ps,  Ts)  +  B((j)  —  (/>),  where  if>s  is  the  mass-specific  solid  phase  free 
energy.  Here,  'ips  is  the  thermal  component  of  free  energy.  Rate-dependent  compaction  is  the 
dominate  dissipation  mechanism  for  compaction  shocks  because  3>  e^. 

The  analysis  indicates  that  increasing  porosity  decreases  the  dissipative  work  rate  within 
shocks  but  increases  the  integrated  dissipative  work  over  shock  rise  times  which  is  indicative 
of  enhanced  sensitivity.  Representative  predictions  for  the  dissipation  rate  and  dissipative 
work  profiles  for  a  pressure  dominated  compaction  shock  corresponding  to  Up  =  500  m/s  are 
shown  in  Fig.  18.  Qualitative  agreement  between  mesoscale  and  continuum  theory  predictions 
is  reasonable,  though  quantitative  discrepancies  exist  due  to  differences  in  their  constitutive 
theories.  Figure  19(a)  illustrates  a  comparison  between  predicted  variations  in  dissipative 
work  with  shock  pressure  and  initial  packing  density;  again,  good  qualitative  agreement  exist 
in  all  cases.  The  relation  between  dissipative  work  and  shock  pressure  can  be  accurately 
described  by  a  density-dependent  power  law.  The  continuum  theory  results  in  shock  thick¬ 
nesses  that  are  comparable  to  those  predicted  by  mesoscale  M&S,  where  the  thickness  is 
largely  independent  of  </>o ,  but  appreciably  decreases  with  increasing  Up,  which  is  consistent 
with  measured  data  for  low  density  HMX.  [12]  A  comparison  between  predicted  shock  rise 
times  is  given  in  Fig.  19(b). 

2.  Rapid  and  successive  inert  shock  loading  of  low  density  HMX  was  computation¬ 
ally  examined.  [Rao  and  Gonthier,  2015  A]  Shock  compaction  of  granular  explosives  can 
trigger  combustion  that  results  in  detonation,  even  for  relatively  mild  shocks.  It  is  well  es¬ 
tablished  that  a  primary  (or  lead)  shock  can  desensitize  the  material  to  subsequent  shocks  by 
reducing  porosity.  This  phenomenon,  referred  to  as  shock  desensitization,  has  been  observed 
during  Deflagration-to-Detonation  Transition  (DDT)  of  low  density  granular  explosives  in 
which  complex  interactions  between  impact  and  combustion-supported  shocks  influence  the 
initiation  process.  In  this  effort,  a  computational  analysis  was  performed  to  characterize  how 
rapid  successive  planar  shock  loading  of  low  density  HMX  affects  dissipation.  Mesoscale  M&S 
was  used  to  predict  effective  shock  profiles  and  to  examine  hot-spot  fields  induced  by  pore 
collapse.  Resolved  shock  profiles  were  compared  to  those  given  by  the  continuum  compaction 
theory,  and  both  were  analyzed  in  a  thermodynamic  space  that  highlights  desensitization 
effects. 

Figures  20(a)  and  (b)  show  the  spatial  variation  in  dissipated  work  iVd  at  various  times  fol¬ 
lowing  initial  impact  predicted  by  mesoscale  M&S  and  the  continuum  compaction  theory, 
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(a) 


(b) 


Figure  18:  A  comparison  between  mesoscale  and  continuum  theory  predictions  for  shock  dissipation 
rate  and  dissipative  work  profiles  corresponding  to  Up  =  500  m/s.  Shock  profiles  are  shown  for 
different  values  of  initial  solid  volume  fraction  0O.  The  shock  is  propagating  from  left  to  right  in 
the  figure. 


(a) 


(b) 


Figure  19:  A  comparison  between  mesoscale  and  continuum  theory  predictions  for  (a)  the  variation 
in  effective  dissipative  work  with  shock  pressure  and  (b)  the  variation  in  compaction  shock  rise 
time  with  supporting  piston  speed. 
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(c) 
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Figure  20:  A  comparison  of  mesoscale  M&S  and  continuum  compaction  theory  predictions  for 
dissipative  work  profiles  induced  by  successive  planar  shock  loading  of  low  density  HMX  having 
4>o  =  0.68.  The  successive  supporting  piston  speeds  are  Up\  =  300  m/s  and  Up2  =  500  m/s. 
(a)  Mesoscale  and  (b)  continuum  theory  spatial  fields  at  various  times.  Fields  are  plotted  in  a 
right  propagating  piston-attached  reference  frame;  the  piston  surface  is  located  at  £  =  0  mm.  (c) 
Mesoscale  and  (d)  continuum  theory  predictions  plotted  in  pressure-dissipated  work  phase  space. 
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respectively,  for  <fi o  =  0.68.  Corresponding  times  differ  between  the  simulations  due  to  dif¬ 
ferences  in  domain  lengths  and  the  time  at  which  the  piston  speed  transitions  from  Up\  to 
Up-2-  In  both  cases,  the  lead  compaction  shock,  which  eliminates  much  of  the  initial  porosity, 
propagates  a  relatively  long  distance  into  the  domain  before  the  piston  speed  increases  to 
U p2  ■  Subsequently,  a  quasi-steady  secondary  shock  quickly  develops  that  eliminates  residual 
porosity  in  material  precompacted  by  the  lead  shock.  The  secondary  shock  is  overdriven 
relative  to  the  ambient  material  which  has  lower  impedance  than  the  denser  material  through 
which  it  propagates.  Consequently,  it  experiences  a  gradual  decay  in  pressure  as  it  over¬ 
takes  the  lead  shock  and  transits  to  the  ambient  material.  The  pressure  decay  predicted  by 
mesoscale  M&S  is  faster  than  that  of  the  continuum  theory  due  to  its  shorter  domain  length 
that  facilitates  rapid  reverberation  of  release  waves  between  the  shock  and  piston.  The  shock 
transmission  process  is  accompanied  by  a  relatively  modest  increase  in  solid  particle  velocity 
and  a  substantial  increase  (~  350%)  in  dissipative  compaction  work  which  mesoscale  sim¬ 
ulations  indicate  can  result  in  ubiquitous  formation  of  hot-spots.  Mesoscale  M&S  and  the 
continuum  theory  predict  qualitatively  similar  features  and  exhibit  reasonable  quantitative 
agreement. 

Figures  20(c)  and  (d)  display  the  predictions  of  (a)  and  (b)  in  a  phase  space  characterized 
by  solid  pressure  Ps  and  Wd ■  This  space  is  chosen  because  Ps  represents  a  measure  of  shock 
strength  and  Wd  represents  a  measure  of  hot-spot  formation  within  the  material.  For  compar¬ 
ative  purposes,  single  shock  end  state  predictions  given  by  both  descriptions  are  also  shown 
in  the  figure.  Mesoscale  M&S  of  single  shocks  indicate  that  appreciable  hot-spot  agglomera¬ 
tion  results  for  approximately  Wd  >  100  kJ/kg  and  Ps  >  1.6  GPa  for  </> o  =  0.68;  this  region 
is  also  highlighted  in  the  figure  for  reference.  Though  speculative  in  the  absence  of  reac¬ 
tive  mesoscale  data,  it  is  reasonable  to  expect  that  shocked  material  having  tightly  packed 
agglomerated  hot-spots  may  significantly  reduce  reactive  hot-spot  interaction  times,  possi¬ 
bly  resulting  in  prompt  ignition  and  initiation.  This  assertion  is  consistent  with  the  almost 
“discontinuous”  initiation  events  observed  by  McAfee,  et  ah, [10]  during  DDT  in  low  density 
HMX  resulting  from  the  interaction  of  piston  and  combustion  driven  compaction  shocks.  Such 
observations  reinforce  the  need  to  better  characterize  dissipation  and  ignition  by  rapid  and 
successive  shock  loading  of  low  density  explosives. 

Identified  on  the  autonomous  solution  trajectories  of  Figs.  (20)  (c)  and  (d)  are  the  equilib¬ 
rium  end  states  behind  the  lead  and  secondary  quasi-steady  shocks,  where  the  lead  shock  end 
states  lie  on  the  single  shock  equilibrium  curves.  Unlike  these  autonomous  solution  trajecto¬ 
ries,  those  associated  with  the  interaction  process  between  the  lead  and  secondary  shocks  are 
time-dependent;  only  a  single  representative  trajectory  is  shown  for  each  of  these  processes 
in  the  plots.  Again,  mesoscale  M&S  and  the  continuum  theory  give  comparable  results.  A 
few  features  are  noteworthy.  First,  the  lead  shock  strength  is  sufficient  to  eliminate  most 
porosity  resulting  in  relatively  small  additional  increases  in  Wd  by  the  secondary  shock.  Con¬ 
sequently,  in  all  cases  the  value  of  Wd  induced  by  the  secondary  shock  is  substantially  lower 
than  that  induced  by  a  single  shock  of  equivalent  pressure  which  is  indicative  of  significant 
shock  desensitization.  Second,  based  on  an  analysis  of  hot-spots,  appreciable  agglomeration 
only  results  as  the  secondary  shock  overtakes  the  lead  shock  for  the  values  of  Up\  and  Up2 
imposed  here.  This  prediction  demonstrates  how  the  strength  of  the  secondary  shock  needed 
to  induce  agglomeration  is  dependent  on  the  strength  of  the  lead  shock  and  the  value  of  </>o- 
Third,  increasing  <po  tends  to  suppress  dissipation  and  enhance  desensitization. 
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3.  A  multiphase  theory  for  the  detonation  of  low  density  explosives  containing  an 
arbitrary  number  of  solid  components  was  formulated.  [Crochet,  et  al.,  2012;  Cro¬ 
chet,  2013;  Crochet  and  Gonthier,  2015]  Multiphase  continuum  models  are  commonly  used  to 
predict  the  shock,  combustion  and  detonation  of  granular  energetic  mixtures  containing  solid 
reactants  and  gaseous  products.  These  models  often  include  phase  interaction  terms  that 
formally  satisfy  the  strong  form  of  the  Second  Law  of  Thermodynamics  and  provide  flexibil¬ 
ity  in  distributing  dissipation  between  phases  arising  from  non-equilibrium  phenomena.  This 
effort  formulated  a  thermodynamically  compatible  constitutive  theory  for  reactive  systems 
containing  an  arbitrary  number  of  solid  components.  The  theory  represents  a  rigorous  ex¬ 
tension  of  the  two-phase  theory  formulated  by  Bdzil,  et  al.  [4],  based  on  the  well-studied 
Baer-Nunziato  model.  Forms  of  the  gas-solid  and  solid-solid  interphase  sources  suggested  by 
general  reactions  of  type  A  — >•  B  are  considered,  where  the  combustion  processes  discussed 
in  [4]  are  treated  as  a  special  case.  The  model  energetics  are  augmented  by  supplemental 
evolutionary  equations  that  track  local  changes  in  phase  temperatures  due  to  dissipative  and 
transport  processes  allowing  for  the  identification  of  dominant  energetic  processes.  This  ca¬ 
pability  provides  a  means  to  identify  system  parameters  (e.g.,  metal  particle  size  and  mass 
fraction  in  metalized  energetic  mixtures)  which  optimize  performance  metrics.  Predictions 
were  obtained  for  mixtures  of  granular  HMX  and  aluminum  to  demonstrate  model  features 
and  to  highlight  the  effect  of  aluminum  particle  self-heating  by  oxidation  on  detonation.  Pre¬ 
dicted  spatial  profiles  for  mechanical  fields,  and  the  heating  contributions  from  individual 
dissipative  processes,  illustrated  how  aluminum  particle  size  can  affect  the  the  coupling  of 
oxidative  heating  to  the  explosive  reaction  zone.  This  extended  theory  is  currently  being 
used  by  the  PI  to  computationally  characterize  shock  initiation  of  HMX-A1  composites. 

4.  A  simple  reactive  burn  model  was  formulated  that  describes  shock  induced  com¬ 
bustion  of  low  density  HMX.  [Rao  and  Gonthier,  2014  B]  The  burn  model  incorporated 
the  ignition  criterion  based  on  effective  hot-spot  formation  rates  discussed  in  Section  3  (Item 
4).  The  effective  rate  was  related  to  bulk  wave  pressure  through  a  power  law.  The  burn  model 
was  implemented  in  the  two-phase  (explosive  and  product  gas)  non-equilibrium  continuum 
theory  given  by  Eq.  (4)  to  computationally  examine  how  ignition  and  initiation  vary  with 
initial  packing  density  and  shock  strength. 

Of  particular  importance  in  this  effort  is  the  volumetric  mass  exchange  term  Cs  that  governs 
the  explosive  combustion  rate.  It  is  common  to  take  Cs  =  Cs(ps,  <f>s,  Pg),  where  the  functional 
dependency  on  gas  pressure  reflects  strand  burn  data.  Because  ignition  results  from  dissipa¬ 
tive  heating  at  the  pore  scale  during  shock  compaction,  it  is  plausible  to  use  an  expanded 
expression  that  accounts  for  the  role  of  hot-spots  in  triggering  ignition. 

The  ignition  induction  period  behind  shocks  depends  on  both  mesostructure  and  shock 
strength.  To  this  end,  Eq.  (4)  is  supplemented  by  a  history  variable  rate  equation  that 
describes  the  evolution  of  a  dimensionless  ignition  variable  A j: 

dh  =  i  (ps-p0y 

dts  tr  \  Pr  /  ’ 

where  tr  and  Pr  are  a  reference  time  and  pressure,  rn  is  a  history  variable  exponent,  and 
Pq  (ps  0  for  this  effort)  is  the  ambient  material  pressure.  The  ignition  variable  is  defined 
for  0  <  Xj  <  1,  where  A/  =  1  denotes  ignition.  The  solid  pressure  Ps  is  chosen  as  the 
history  variable  because  it  drives  compaction  induced  dissipation  which  represents  a  primary 
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hot-spot  formation  mechanism  in  the  continuum  theory.  Assuming  a  constant  shock  pressure 
and  taking  ts  =  0  at  the  shock  front,  Eq.  (11)  can  be  directly  integrated  to  give  A /  = 
(' ts/tr)((Ps  —  Po)/Pr)m ,  where  ts  is  time  since  shock  passage  for  the  explosive  phase.  For 
ignition,  A/  =  1  and  ts  =  r*;  also,  Ps  «  P  because  (j)s  «  1  in  the  compacted  state  for  shocks 
considered  in  this  effort.  Thus,  using  the  criterion  P"t*  =  ^  from  Section  3  (Item  4),  and 
taking  m  =  n,  the  expression  1  =  (r* /tr)((Ps  —  Po)/Pr)n  is  obtained  enabling  Eq.  (11)  to  be 
expressed  as 

d\ ,L=[  fr( Ps-PoT  0  <  A,  <  1 

1  0  otherwise  ’ 

where  the  constant  parameters  a,  I*,  and  n  are  mesostructure-dependent.  The  volumetric 
mass  exchange  term  is  then  taken  to  be 


0  for  0  <  A/  <  1 

apscj)sPg  for  A/  =  1, 

where  the  constant  a  is  a  burn  prefactor.  Here,  it  is  assumed  that  negligible  gas  is  pro¬ 
duced  during  the  induction  period  corresponding  to  0  <  A/  <  1  as  hot-spots  cook-off  within 
the  microstructure.  It  is  further  assumed  that  the  induction  period  is  energetically  incon¬ 
sequential  at  the  macro-scale.  This  term  assumes  that  vigorous  combustion  progresses  at  a 
rate  determined  by  the  gas  pressure  at  the  instant  of  ignition;  consequently,  this  model  can 
give  predictions  that  differ  from  reactive  burn  models  that  impose  a  constant  ignition  pres¬ 
sure.  Different  modes  of  transition  to  detonation  are  not  expected  to  be  accurately  described 
by  Eq.  (13)  due  to  the  simplistic  nature  of  the  pressure-dependent  burn  law.  This  simple 
model  is  chosen  as  a  first  attempt  to  better  understand  how  the  ignition  period  can  influence 
detonation  transition  over  a  wide  range  of  pressure. 


Figure  21  summarizes  predictions  for  the  variation  in  time  t*  and  distance  x*  to  detonation 
break-out  within  the  material  following  impact  with  input  shock  pressure  for  porous  HMX 
having  (f>0  =  0.68  and  0.84;  these  plots  are  commonly  referred  to  as  Pop-plots.  Also  shown  in 
the  figure  are  initiation  time  and  run  distance  data  reported  in  the  literature  for  porous  HMX 
having  (j)0  =  0.65 [7,  8,  10,  12]  and  (f>0  =  0.86.  [8]  Though  these  materials  have  similar  packing 
densities,  the  particle  size  distributions  differ;  therefore,  these  results  are  only  intended  to 
illustrate  the  leading-order  effects  of  packing  density  on  ignition  and  initiation.  As  seen 
in  Fig.  21(a),  the  initiation  time  for  each  material  (dotted  curve)  is  largely  controlled  by 
ignition  for  strength  dominated  shocks,  where  ignition  time  curves  are  plotted  as  solid  lines; 
such  processes  are  therefore  ignition  controlled.  These  ignition  curves  are  slightly  offset  from 
those  predicted  by  the  continuum  theory  at  higher  pressure  (dashed  lines)  due  to  the  influence 
of  shock  rise  time  in  the  implementation  of  Eq.  (12).  As  input  pressure  increases,  predicted 
initiation  times  are  long  compared  to  ignition  times  and  are  controlled  by  the  growth  of 
reaction  to  detonation;  such  processes  are  growth  controlled.  Predicted  initiation  times  for 
(j>0  =  0.68  compare  favorably  to  those  measured  by  Dick[7]  for  high  pressure  shocks.  The 
shock  pressure  required  for  transition  from  ignition  to  growth  controlled  initiation  increases 
with  packing  density,  and  is  estimated  to  be  approximately  0.2  GPa  for  (j)0  =  0.68  and  0.4  GPa 
for  (f) o  =  0.84. 


Figure  21(b)  summarizes  the  corresponding  predicted  variation  in  run  distance  to  detonation. 
A  couple  observations  are  noteworthy.  First,  due  to  a  complex  transition  mode,  ignition  con¬ 
trolled  initiation  is  characterized  by  run  distances  associated  with  the  evolution  of  detonation 
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Figure  21:  Plots  showing  the  (a)  ignition  and  initiation  times,  and  (b)  initiation  distances  for 
<f>0  =  0.68  and  <p0  =  0.83. 


in  compacted  material  behind  the  lead  shock  and  the  break-out  of  detonation  into  ambient 
material.  These  distances  merge  with  increasing  shock  pressure  as  the  initiation  process  be¬ 
comes  growth  controlled.  Second,  slopes  of  the  initiation  curves  (given  by  Ologa^/cHog-P) 
slightly  decrease  with  increasing  shock  pressure  for  ignition  controlled  processes,  but  approach 
an  approximately  constant  value  for  growth  controlled  processes  that  depends  on  0O.  Higher 
<fi0  results  in  a  marginally  higher  slope,  possibly  due  to  the  larger  change  in  effective  hot-spot 
formation  rate  with  pressure.  Also  shown  in  Fig.  21(b)  for  perspective  is  the  measured  run 
distance  to  detonation  for  the  plastic-bonded  explosive  PBX  9404  (98%  TMD  HMX).  Relative 
to  <f> o  =  0.84,  the  run  distance  curve  for  this  material  has  higher  slope,  suggesting  a  larger 
variation  in  hot-spot  formation  rate  with  pressure,  and  requires  higher  shock  pressure  for  a 
given  run  distance;  as  such,  it  is  considered  less  sensitive  than  low  density  HMX. 

Though  the  reactive  burn  model  formulated  in  this  effort  gives  results  that  are  commensurate 
with  time  and  run  distance  data,  not  all  details  of  the  transition  mechanism  are  accurately 
predicted  following  ignition  due  to  the  simple  form  of  the  pressure-dependent  burn  law  used. 
This  deficiency  is  particularly  apparent  for  ignition  controlled  initiation  that  involves  complex 
temporal  and  spatial  interactions  between  compaction  and  combustion  waves  within  material 
behind  the  lead  shock. 

5.  Computational  analysis  of  weak  initiation  of  low  density  HMX.  [Rao  and  Gonthier, 
2015  B;  Rao  and  Gonthier,  2016  A;  Rao  and  Gonthier,  2016  B]  Low  pressure  (weak)  initi¬ 
ation  of  low  density  granular  HMX  occurs  by  a  complex  mechanism  that  leads  to  a  prompt 
(or  discontinuous)  transition  to  detonation  within  the  material  due  to  compaction  shock  in¬ 
teractions.  These  interactions  influence  ignition,  flame  spread,  and  subsequent  transition 
by  affecting  dissipative  heating  within  the  microstructure  during  pore  collapse.  Details  of 
the  transition  mechanism  depend  on  initial  packing  density  and  input  shock  strength.  In 
this  effort,  computations  were  performed  using  a  multiphase  reactive  flow  model  to  exam¬ 
ine  how  the  transition  mechanism  varies  with  input  shock  strength  for  granular  HMX  (68% 
TMD).  As  mentioned,  the  burn  model  outlined  in  Item  4  of  this  section  does  not  accu- 
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rately  describe  details  of  the  transition  mechanism  due  to  the  overly  simplistic  nature  of 
the  pressure-dependent  burn  law  despite  its  accuracy  in  predicting  Pop-plots.  The  improved 
model  accounts  for  pressure-dependent  ignition,  and  subsequent  burn  depends  on  the  local 
dissipative  work,  porosity,  and  pressure.  The  dependence  on  dissipative  work  is  motivated 
by  mesoscale  M&S  that  indicated  a  significant  increase  in  hot-spot  size  and  spatial  proximity 
within  the  microstructure  as  shock  induced  dissipative  work  increases,  suggesting  an  increase 
in  flame  spread  rate.  Predictions  highlight  the  variation  in  transition  mechanism  with  in¬ 
creasing  input  shock  strength  and  identify  conditions  that  favor  the  formation  of  retonation 
waves  during  transition. 


The  volumetric  mass  exchange  term  Cs  is  again  described  by  Eq.  (13),  where  the  subsequent 
burn  associated  with  gas  production  and  energy  release  is  now  given  by 


aL  (  do  )  <Tre/  (  Prlf  )  f°r  (t)s  —  (f>crit 

((XL  +  aH)  ()|j)  ef  (7^7)  for  >  4>crit ■ 


(14) 


Here,  a  represents  a  dimensionless  volumetric  burn  surface  area  prefactor,  do  is  particle  size, 
and  aref  is  a  characteristic  regression  velocity  corresponding  to  pressure  Pref-  Generally, 
a  =  a(4>s,Wd),  where  is  mass  specific  shock  dissipative  work  that  can  account  for  desen¬ 
sitization.  This  burn  rate  reflects  a  transition  from  normal  regressive  burn,  characterized  by 
aL,  to  hot-spot  facilitated  burn,  characterized  by  an,  for  strong  pressure  dominated  shocks 
that  have  (f)s  ~  1  and  high  values  of  wd-  Hot-spot  facilitated  burn  is  associated  with  a  signif¬ 
icant  increase  in  burn  area  due  to  the  ubiquitous  formation  of  hot-spots.  The  dependence  on 
(ficrit  is  chosen  because  DDT  tube  tests  [10]  indicate  a  drop  in  burn  rate  for  approximately 
(j>s  >  0.95  prior  to  the  onset  of  a  hot-spot  facilitated  burn-type  mechanism.  For  (j)s  <  4>crit,  a 
constant  prefactor  is  used  and  the  burn  is  solely  dependent  on  gas  pressure.  The  prefactor 
an  is  given  by 


an  =  ar 
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In  this  effort,  the  values  are  do  =  100  /irn,  aref  =  10  m/s,  Pref  =  0.1  GPa,  aL  =  0.001, 
amin  =  0.001,  amax  =  0.2,  wref  =  20  kJ/kg,  and  A  =  0.5  kJ/kg.  Different  detonation 
transition  mechanisms  can  be  described  by  Eq.  (14)  due  its  dependence  on  4>s,  Pg,  and  wd- 
Though  conceptually  similar  to  other  models  [13,  14],  this  model  incorporates  microstructure- 
dependent  hot-spot  predictions  obtained  by  meso-scale  M&S  and  accounts  for  hot-spot  facili¬ 
tated  burn  in  terms  of  compaction  shock  dissipation  to  better  characterize  desensitization  by 
shock  interaction  effects. 


Figure  22  illustrates  a  predicted  transition  mechanism  in  low  density  HMX  having  <p 0  =  0.68 
for  an  input  shock  corresponding  to  a  piston  speed  of  Up  =  200  m/s.  The  input  shock  in 
this  case  is  strength  dominated  and  it  triggers  an  ignition  controlled  initiation  process.  The 
steady  primary  (input)  compaction  shock,  labeled  ci,  eliminates  most  of  the  porosity  as  it 
propagates  a  relatively  long  distance  into  the  ambient  material  prior  to  the  onset  of  burn 
near  the  piston  surface  at  approximately  t  =  25  /is  following  impact.  Subsequently,  regressive 
burn  near  the  piston  drives  the  formation  of  a  secondary  compaction  shock,  labeled  C2, 
due  to  the  low  permeability  of  combustion  gases  into  dense  material  precompacted  by  the 
primary  shock.  The  secondary  shock  eliminates  residual  porosity  resulting  in  the  formation 
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<t>Q  =  0.68 
Up  =  200  m/s 

Cj  -  Primary  Compaction  Wave 
c2  —  Secondary  Compaction  Wave 
b  -  Ignition  Bum  Front 
D  -Detonation 
R  -Retonation 
S  -  Shock 
Sj- Inert  Shock 


t  (|4S) 


(a) 


(b) 


Figure  22:  Predicted  (a)  bulk/mixture  pressure  contours  showing  compaction  shock  and  ignition 
front  trajectories  in  the  characteristic  plane;  (b)  solid/explosive  volume  fraction  contours  in  the 
characteristic  plane;  (c)  variation  in  bulk  pressure  with  position  at  various  times;  and  (d)  variation 
in  dissipative  compaction  work  with  position  at  various  times.  The  material  is  low  density  HMX 
(cj) o  =  0.68)  and  the  input  shock  corresponds  to  a  supporting  piston  speed  of  Up  =  200  m/s. 
Position  is  shown  in  a  piston-attached  reference  frame. 
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of  an  inert  solid  “plug”  that  increases  in  width  as  it  propagates  forward  [10].  This  burn 
supported  secondary  shock  causes  less  dissipation  than  those  induced  by  weaker  input  shocks 
(~  Up  <  150  m/s)  because  of  slightly  lower  residual  porosity  behind  the  primary  shock 
which  is  indicative  of  desensitization.  The  secondary  shock  overtakes  the  primary  shock  at 
approximately  t  =  35  /xs  causing  a  significant  increase  in  dissipation  with  a  slight  reduction 
in  pressure.  Consequently,  the  rapid  onset  of  vigorous  burn  is  predicted  to  occur  within  the 
compacted  region  at  approximately  t  =  43  /is  after  a  short  ignition  delay  following  shock 
passage.  At  this  location  and  time,  Wd  — >•  wref  =  20  kJ/kg  and  an  quickly  increases  resulting 
in  hot-spot  facilitated  burn.  This  rapid  onset  of  vigorous  burn  results  in  a  “thermal-explosion¬ 
like”  (or  SDT-like)  event  which  produces  a  forward  propagating  reactive  shock  and  a  backward 
propagating  inert  shock.  The  reactive  shock  quickly  transitions  to  a  Chapman- Jouguet  (CJ) 
detonation,  whereas  the  inert  shock  slowly  increases  in  strength  until  it  abruptly  transitions 
into  a  strong  reactive  compression/retonation  wave  as  Wd  — >  wref.  These  predictions  indicate 
that  the  formation  of  retonation  waves  requires  the  occurrence  of  a  SDT-like  event  which 
drives  a  backward  propagating  inert  shock  that  becomes  sufficiently  strong  to  trigger  hot¬ 
spot  facilitated  burn  before  interacting  with  the  primary  burn  front. 

Though  not  shown  here,  weaker  input  shocks,  which  also  involve  SDT-like  events,  do  not  result 
in  retonation  as  the  inert  shock  strength  is  insufficient  to  trigger  hot-spot  facilitated  burn  prior 
to  interacting  with  the  primary  burn  front.  Stronger  input  shocks  that  are  pressure  dominated 
eliminate  all  porosity  and  result  in  a  continuously  accelerating  combustion/detonation  wave 
that  forms  within  the  compacted  region  whose  pressure  exceeds  the  Chapman-Jouguet  value 
of  the  ambient  material.  In  these  cases,  the  transition  mechanism  is  largely  continuous,  with 
ignition  induction  time  comparable  to  that  required  for  transition  following  ignition. 
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